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A  Technique  to  Assess  the  Characteristics  of  Bottom  and  Subbottom  Marine 
Sediments  (TR  DRP-95-3) 


ISSUE:  The  accurate  characterization  of  bot¬ 
tom  and  subbottom  sediments  is  critical  in  all 
aspects  of  dredging  operations  from  planning 
to  actual  removal  and  disposal  of  the  s^iment. 

A  need  existed  to  develop  a  rapid  geophysical 
technique  for  determination  of  the  physical 
properties  of  materials  scheduled  for  dredging. 

RESEARCH:  A  decision  was  made  that 
acoustic  principles  might  be  the  best  to  use  in 
developing  a  reliable  technique  to  determine 
the  dredgeability  of  sediment.  The  work  had 
two  objectives: 

•  Develop  an  electronic  package  to  send  and 
analyze  acoustic  signals  to  provide  geo¬ 
physical  information  from  the  acoustic  re¬ 
flectivity  strength  of  the  signals. 

•  Survey  existing  geophysical  equipment 
and  develop  a  rapid  technique  and  appro¬ 
priate  instrumentation  to  determine  the 
physical  properties  of  bottom  and  sub¬ 
bottom  materials. 


n.  Continued  field  testing  at  a  second  site. 

in.  Incorporation  of  results  of  reimbursable 
field  tests  and  continued  development 
of  acoustic  impedance  (AI)  approach. 

SUMMARY:  The  study  resulted  in  an  optimum 
AI  system.  The  report  includes  discussions  of 
equipment  and  technical  development  of  the  AI 
concept;  details  of  survey  planning;  and  data  pro¬ 
cessing,  interpretation,  and  presentation.  Since 
several  AI  surveys  were  conducted  as  reimburs¬ 
able  projects  during  the  course  of  the  research 
study,  selected  sites  were  included  in  the  report 
to  illustrate  the  AI  concept,  its  application,  and 
ground-tmth  comparisons. 

AVAILABILITY  OF  REPORT:  The  report  is 
available  on  Interlibrary  Loan  Service  from  the 
U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  Library,  telephone  (601)  634- 
2355.  National  Technical  Information  Service 
(NTIS)  report  numbers  may  also  be  requested 
from  the  WES  librarians. 


The  effort  was  accomplished  in  three  separate 
but  consecutive  phases: 

I.  A  field  demonstration  of  an  acoustic 
profiler’s  ability  to  measure  and  relate 
acoustic  reflectivity  to  density. 


To  purchase  a  copy,  call  NTIS  at  (703)  487- 


4650. 
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Summary 


This  report  presents  the  theoretical  concept,  assembly,  and  field  testing  of  a 
waterborne  seismic  acoustic  impedance  technique  which  has  been  developed  to 
characterize  bottom  and  subbottom  sediments  as  they  relate  to  removal  by 
dredging.  This  method,  developed  under  the  Dredging  Research  Program 
(DRP),  provides  estimates  of  in  situ  density  and  soil  type  in  a  rapid,  cost- 
effective  manner  using  digital  acoustic  subbottom  profiling  methodology.  In 
situ  densities  obtained  by  the  acoustic  impedance  technique  to  date,  when  com¬ 
pared  to  those  obtained  by  conventional  means  at  several  different  sites  under 
a  wide  variety  of  marine  conditions,  have  statistically  been  within  ±  10  per¬ 
cent.  However,  only  marine  sediments  considered  to  be  fully  saturated,  inor¬ 
ganic,  and  uncontaminated  have  been  investigated  and  “ground  truthed.”  After 
comparisons  with  ground  truth  information  and  laboratory  testing,  a  critical 
analysis  of  the  acoustic  impedance  technique  reveals  it  to  be  a  valid  and  useful 
approach  to  bottom  and  subbottom  material  and  density  prediction.  While 
some  development  is  still  needed  to  fully  establish  advantages  and  limitations, 
its  potential  usefulness  warrants  technology  transfer  now,  provided  proper  cau¬ 
tions  are  observed.  This  supposition  is  corroborated  by  the  fact  that  numerous 
reimbursable  surveys  have  been  successfully  conducted  while  products  were 
stiU  in  the  development  stage.  Each  site  surveyed  provided  valuable  input  to 
the  research  and  development  evolutionary  processes  and  enabled  researchers 
to  fine-tune  procedures  while  still  providing  a  useful  and  timely  service  to 
sponsoring  Districts  and  Divisions. 

The  following  benefits  have  been  derived  as  a  result  of  this  work: 


Primary  Products 

a.  A  rapid  geophysical  technique  using  acoustics  to  characterize 
bottom/subbottom  conditions. 

b.  A  computer  visualization  process  for  displaying  subbottom  density  and 
material  type  predictions  and  volumetric  calculations  in  three 
dimensions. 


c.  Specifications  and  assembly  of  waterborne  integrated  geophysical  system 
(assembly  and  components  funded  by  sources  other  than  DRP). 

Support  Products 

a.  Sediment  density  estimates. 

b.  Material  type  predictions. 

c.  Subbottom  stratigraphic  infonnation  for  project  planning  and  design. 

d.  Project  monitoring. 

e.  Volume  estimates  of  sediment. 

/.  Contractual  evidence. 

g.  Geologic  database  for  long-tenn  planning. 

Areas  of  Potential  Cost  Savings 

(The  following  are  based  upon  DRP  benefits  analysis  (Griffis  1994).) 

a.  Direct  savings: 

(1)  Continuous  coverage  of  bottom  and  subbottom  conditions. 

(2)  Optimum  boring  placement  and  elimination  of  those  unnecessary. 

(3)  Accurate  assessment  of  material  type,  density,  and  volume. 

b.  Environmental  enhancements'. 

(1)  Presence  of  cultural  resources. 

(2)  Long-tenn  monitoring  of  disposal  benns. 

c.  Mission  enhancements: 

(1)  Location  of  under-channel  utilities. 


(2)  Detennination  of  fluid  mud  zone  thickness. 


d.  Cost  avoidances: 


(1)  Contractor  claims  on  volumes  of  material  and  surprises  associated 
with  material  changes. 

(2)  Elimination  of  needless  borings. 

(3)  Change  of  chaimel  routing  to  avoid  undesirable  conditions. 

Spinoffs  and  Future  Prospects 

(The  following  are  authors’  opinions  based  upon  observations  made  at  some 
15  sites  surveyed  and  communication  with  product  users.) 

a.  Spinoffs: 

(1)  Incorporation  of  airbome/waterbome  Navy  EM  for  reconnaissance 
determinations  of  sediments  and  topographic/bathymetric  surveying 
(large  area  site  characterization). 

(2)  Determination  of  navigable  depth  (fluid  mud  assessment). 

(3)  Assessment  of  integrity  of  disposal  bemis. 

b.  Future  prospects: 

(1)  Location  of  contaminated  materials. 

(2)  Location  of  marine  sand  deposits. 

(3)  Location  of  dredging  hazards  and  utility  crossings. 

(4)  Bottom/subbottom  conditions  in  lakes  and  reservoirs. 

(5)  Description  of  geologic  foraiations. 


Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

cubic  yards 

0.7645549 

cubic  meters 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

meters 

inches 

0.0254 

meters 

miles  (U.S.  statue) 

1,609.347 

meters 

tons  (force) 

8,896.444 

newtons 

yards 

0.9144 

meters 

1  Introduction 


Background 

During  fiscal  year  1988  the  U.S.  Amiy  Coips  of  Engineers  (USAGE)  began 
a  research  and  development  program  of  major  proportions  to  address  problems 
and  needs  arising  in  the  performance  of  its  dredging  mission.  This  effort 
ealled  the  Dredging  Research  Program  (DRP),  is  a  major  thrust  toward  devel¬ 
oping  improved  technologies  that  can  reduce  the  costs  of  dredging  operations. 
The  program  concentrates  on  problem  areas  related  to  the  physical  aspects  of 
dredging  or  dredging  projeets. 

The  USAGE  Directorate  of  Civil  Works  is  involved  in  virtually  every  navi¬ 
gation  dredging  operation  perfonned  in  the  United  States.  The  Corps  dredging 
mission  entails  maintenance  and  improvement  of  25,000  miles*  of  commer¬ 
cially  navigable  channels  serving  400  ports,  including  130  of  the  nation’s 
150  largest  cities.  Connecting  waterways  to  the  nation’s  ports  and  harbors 
handle  about  2  billion  tons  of  commerce  each  year  as  waterborne  transport 
continues  to  be  the  most  cost-  and  energy-efficient  means  of  shipping  bulk 
cargos.  Additionally,  the  waterways  network  constitutes  an  infrastructure 
component  essential  to  the  nation’s  defense  capabilities. 

The  Corps  of  Engineers  dredges  an  average  of  250  to  300  million  cu  yd  of 
sedimentary  material  at  a  current  expenditure  level  of  about  $400  million  per 
year  in  order  to  maintain  and  operate  the  nation’s  existing  navigation  system 
Dredging  is  the  single  most  costly  item  in  the  Coips’  Civil  Works  Operations/ 
Maintenance  budget. 


The  Corps  has  been  actively  involved  in  dredging  for  over  160  years. 
Originally,  the  Corps  maintained  a  large  fleet  of  Its  own  dredges,  but  in  recent 
years  the  bulk  of  field  operations  have  shifted  from  the  govermnent  fleet  to  the 
private  sector.  In  order  to  fulfill  its  mission  within  budget  and  time  con¬ 
straints,  the  Corps  will  be  continually  challenged  to  meet  dredging  needs  in  an 
efficient,  cost-effective  manner.  To  meet  these  needs,  it  was  recognized  that  a 
research  and  development  program  would  be  needed  to  identify  areas  witli 


*  A  table  of  factors  for  converting  non-Sl  units  of  measurement  to  SI 
page  XV. 
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high  potential  for  large  paybacks  through  the  investment  of  Corps  research 
dollars.  Eventually,  a  direction  was  established  using  informational  input  from 
those  involved  in  Corps  dredging  operations  which  led  to  the  identification  of 
five  specific  problem  areas.  Accomplishments  reported  herein  fall  within  the 
purview  of  Problem  Area  2,  Material  Properties  Related  to  Navigation  and 
Dredging.  The  specific  objectives  of  Problem  Area  2  are:  (1)  development  of 
instrumentation  and  operating  procedures  for  rapid  surveys  of  fluid  mud  prop¬ 
erties,  (2)  definition  of  navigable  depth  in  fine-grained  sediment,  (3)  develop¬ 
ment  of  instrumentation  for  analyzing  properties  of  consolidated  sediments, 
and  (4)  establishment  of  dredging-related  soil  and  rock  descriptors.  Work 
Unit  32470,  Rapid  Measurements  of  Properties  of  Consolidated  Sediments,  was 
undertaken  to  meet  objective  (3)  cited  above.  To  date,  accomplishments  show 
outstanding  promise,  and  continued  future  development  will  enable  Corps 
personnel  to  more  accurately  assess  the  char^teristics  of  bottom  and  subbot¬ 
tom  sediments  as  they  relate  to  removal  by  dredging. 


Objective 

The  objective  of  DRP  Work  Unit  32470  was  to  develop  a  rapid  geophysical 
technique  and  identify/recommend  the  appropriate  instrumentation  necessary 
for  determination  of  the  physical  properties  of  bottom  and  subbottom  materials 
scheduled  for  dredge  removal. 


Approach 

In  the  early  stages  of  the  DRP,  Work  Unit  32470  was  conceived  to  combat 
a  problem  identified  as  the  Corps’  need  to  categorize  and  quantify  the  dredge- 
ability  of  bottom  sediment  in  a  rapid,  widely  accepted  manner.  At  that  time  it 
was  decided  that  acoustic  principles  might  best  be  used  to  fulfill  the  solution. 
Descriptive  subbottom  information  would  be  valuable  in  dredging  applications 
and  in  preparing  government  estimates  for  dredging  contracts,  selecting  the 
most  efficient  and  productive  dredge,  and  selecting  appropriate  disposal  sites. 
The  full  objective  was  to  develop  an  electronic  package  to  send  and  analyze 
acoustic  signals  to  provide  geophysical  information  such  as  density,  shear 
strength,  and  grain  size  from  the  acoustic  reflectivity  strength  of  the  signals. 
Secondly,  the  program  incorporated  a  survey  of  existing  geophysical  equip¬ 
ment  and  development  of  a  rapid  geophysical  technique  and  the  appropriate 
instrumentation  to  determine  the  physical  properties  of  bottom  and  subbottom 
materials  for  dredge  removal  purposes. 

The  effort  was  divided  into  three  separate  but  consecutive  phases.  Phase  I 
consisted  of  a  field  demonstration  of  an  acoustic  subbottom  profiler’s  ability  to 
measure  and  relate  acoustic  reflectivity  of  bottom  sediments  to  density  of  the 
material.  Phase  II  was  structured  to  involve  a  continuation  of  field  tests  at  a 
second  well-documented  site  while  industry-wide  searches  were  conducted  to 
determine  those  properties  of  sediments  which  were  descriptive  of  the  relative 
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ease  of  removal  by  dredging.  Phase  III  was  originally  to  consist  of  an  evalu¬ 
ation  of  other  seismic  methods,  while  configuring  an  optimum  waterborne  sys¬ 
tem  for  a  variety  of  USAGE  applications.  However,  since  a  high  degree  of 
success  was  achieved  during  Phase  n  and  at  several  follow-on  test  sites  where 
surveys  were  performed  on  a  reimbursable  basis,  it  was  decided  to  incorporate 
representative  data  and  limit  further  evaluations  to  various  types  of  seismic 
sources  while  concentrating  on  full  development  of  the  Acoustic  Impedance 
(AI)  approach.  Further,  through  other  sources,  funds  were  secured  to  purchase 
a  vessel  and  a  full  complement  of  geophysical  equipment  In  the  life  of  this 
program,  the  optimum  AI  system  was  bom  and  rapid  progress  made  in  survey 
performances  and  continued  technological  improvements  in  data  processing 
and  displays.  That  system  is  now  available  for  use  in  a  wide  range  of  problem 
applications,  whether  civil  or  military. 


Scope 

This  report  contains  discussions  on  equipment,  technical  development  of  the 
acoustic  impedance  concept  data  processing  and  interpretation,  data  presenta¬ 
tion  (visualization),  survey  planning,  limitations,  and  recommendations.  Since 
numerous  AI  surveys  have  been  successfully  conducted,  selected  sites  will  be 
used  to  illustrate  the  AI  concept  its  applications,  and  ground  truth  compari¬ 
sons.  This  report  is  not  intended  to  serve  as  a  user’s  guide  or  as  an  Engineer 
Manual.  Procedures  developed  and  reported  herein  require  hands-on  training 
of  specialized  personnel  skilled  in  the  application  of  geophysical  methods,  data 
processing,  and  interpretation. 
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Overview 

The  AI  method  was  developed  to  rapidly  and  accurately  assess  engineering 
properties  of  shallow  marine  sediments,  providing  virtually  continuous  cover¬ 
age  for  delineation  of  both  horizontal  and  vertical  extents  of  those  sediments  as 
they  relate  to  dredging.  Another  goal  was  to  provide  this  “hardware/software 
package”  on  an  affordable,  standardized  computer  platform,  consequently  mak¬ 
ing  the  technology  accessible  to  a  wide  variety  of  users.  Finally  and  most 
importantly,  development  of  “black  boxes”  was  avoided.  As  with  many  geo¬ 
physical  techniques,  there  exists  a  variety  of  situations  or  environments  whose 
response  characteristics  to  a  particular  external  event,  such  as  seismic  excita¬ 
tion,  are  neither  theoretically  nor  empirically  uniquely  defined.  Without  proper 
physical  verification,  or  ground  truthing,  significant  errors  in  interpretation  are 
possible.  In  other  words,  the  state-of-the-art  has  not  advanced  to  the  point 
where  any  one  device  or  system  is  capable  of  accurately  assessing  all  feasible 
situations  because  certain  combinations  of  physical  factors  can  cause  nearly 
identical  seismic  signatures. 

The  AI  method  incorporates  specially  designed  software  packages,  hardware 
configurations,  and  standardized  field  data  collection  methods  into  a  structured 
geoacoustic  modeling  program.  A  geoacoustic  model  as  defined  by  Hamilton 
(1980)  is  “a  model  of  the  real  seafloor  with  emphasis  on  measured,  extrapo¬ 
lated,  and  predicted  values  of  those  properties  important  in  underwater  acous¬ 
tics  and  those  aspects  of  geophysics  involving  sound  transmission.”  Hamilton 
further  shows  that  a  geoacoustic  model  should  include  the  following  informa¬ 
tion  for  each  sediment  layer; 

a.  Sediment  information  (cores). 

b.  Thickness  of  layers. 

c.  Properties  of  overlying  water  mass. 

d.  Locations,  thickness,  and  properties  of  reflectors  within  the  sediment 
body  at  various  frequencies. 
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e.  Properties  of  rock  layers. 

/.  Details  of  the  bottom  topography. 

It  is  this  concept  of  geoacoustic  modeling  that  was  used  to  direct  the  devel¬ 
opment  of  the  AI  method,  specifically  targeting  the  needs  of  the  DRP.  For 
example,  any  particular  application  of  the  AI  method  is  basically  divided  into 
four  phases: 

a.  Data  collection. 

b.  Data  processing. 

c.  Direct  sampling. 

d.  Geoacoustic  modeling. 

Figure  1  is  a  typical  project  flowchart  from  a  comprehensive  AI  study  con¬ 
ducted  off  the  coast  of  Delaware  for  the  Philadelphia  district.  Discussion  of 
the  AI  method  will  begin  with  a  presentation  of  technical  issues  regarding  the 
capability  of  performing  quantitative  analysis  on  acoustic  reflection  data  for  the 
purpose  of  characterizing  engineering  properties  of  marine  sediments.  The 
remaining  discussion  will  follow  according  to  the  above  phases.  Specific  the¬ 
ories,  hardware,  software,  or  procedures  will  be  discussed  in  detail  under  the 
appropriate  topic. 


Acoustic  Impedance  Concept 

The  AI  method  is  an  extension  of  techniques  developed  by  Caulfield  and 
Yim  (1983)  and  Caulfield,  Caulfield,  and  Yim  (1985).  It  must  be  emphasized 
that  the  model  is  an  empirical  technique  which  compensates  for  absorption  in 
each  layer  as  a  function  of  the  center  frequency  of  a  band-limited  seismic 
trace,  corrects  for  spherical  spreading,  and  utilizes  classical  multilayer  reflec¬ 
tive  mathematics  to  compute  reflection  coefficients  at  sediment  horizons. 
Reflection  coefficients  are  converted  to  impedances  and  classified  according  to 
established  relationships  between  acoustic  impedance  and  geotechnical  proper¬ 
ties  of  marine  sediments,  thereby  classifying  the  lithostratigraphy.  Figure  2 
illustrates  the  general  processing  steps  required  by  the  method  in  practice. 
Figure  3  shows  a  hypothetical  echo  sounding  or  seismic  pulse  impinging  on 
the  bottom  and  subbottom.  The  pulse  has  a  finite  time  width  (At)  and  will 
have  different  return  amplitudes  based  on  the  acoustic  impedance  (Z)  and  the 
absorption  (a)  of  each  layer. 


Reference  units  and  the  decibel  scale 

Since  nearly  all  computations  and  results  are  given  in  decibels  (abbreviated 
as  db),  a  brief  explanation  of  the  system  is  provided.  The  unit  of  intensity  in 
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underwater  sound  for  purposes  of  this  report  is  the  intensity  of  a  plane  wave 
having  an  rms  pressure  equal  to  1  dyne  per  square  centimeter  (dyne/cm^ ). 

This  differs  from  the  new  American  National  Standard^  of  1  micropascal 
(1  pPa  =  10'^  dyne/cm^ ).  The  dyne/cm^  standard  is  used  herein  to  correlate 
with  the  work  of  previous  researchers. 

Historically,  decibel  units  have  been  used  for  reckoning  acoustic  quantities 
due  to  their  convenience  in  dealing  with  large  changes  in  variables,  thereby 
simplifying  the  computational  process.  TTie  decibel  is  a  logarithmic  expression 
of  the  ratio  of  two  amounts  of  power,  specifically 

/I  =  10  logioiP/P^) 


where 

P  =  amount  of  power  being  considered 

Pg  =  reference  power 

This  equation  is  used  whenever  a  unit  of  power  or  total  energy  is  being 
considered,  such  as  the  source  level  of  an  acoustic  device  which  is  normally 
expressed  in  kilowatts,  or  when  evaluating  the  total  energy  of  a  wavelet  in  the 
frequency  domain.  In  the  case  of  acoustic  sound  pressures,  the  decibel  value 
is  expressed  as 

«  =  20  log^o(p/Po) 


where 

p  =  given  pressure  level  (db) 

=  reference  pressure  (db) 

It  should  be  emphasized  that  the  decibel  is  a  comparison  of  intensities 
rather  than  a  direct  comparison  of  acoustic  pressures.  The  statement  “20  db 
referenced  to  one  dyne/cm^”  actually  means  “20  db  relative  to  the  intensity  of 
a  plane  wave  of  pressure  equal  to  one  dynetcn?.”  Nearly  all  computations 
involved  in  the  AI  method  are  accomplished  using  the  decibel  scale. 
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Principles  of  Acoustic  Reflection 


The  principles  governing  acoustic  reflection  are  well-known,  are  presented 
in  many  texts,  and  will  not  be  repeated  here.  The  reader  is  referred  to  the 
work  of  Officer  (1958)  for  an  excellent  discussion  of  the  theories  of  sound 
transmission.  Because  AI  is  a  quantitative  approach,  a  discussion  concerning 
the  basic  theory  utilized  in  remotely  determining  the  acoustic  impedance  of  a 
sediment  layer  using  sound  is  appropriate. 

Consider  wave  propagation  across  a  boundary.  If  a  seismic  wave  propagat¬ 
ing  through  a  medium  arrives  at  the  boundary  of  another  medium,  part  of  the 
energy  of  the  wave  will  be  reflected,  part  transmitted,  and  a  portion  absorbed 
within  the  upper  medium.  Starting  from  Snell’s  Law  governing  refraction  and 
applying  the  appropriate  boundary  conditions  of  continuity  of  nonnal  and  tan¬ 
gential  stresses,  it  is  possible  to  represent  the  relationship  between  the  incident, 
reflected,  and  transmitted  waves  in  tenns  of  their  velocities  and  angle  of  propa¬ 
gation  relative  to  the  nonnal.  This  relationship  is  shown  by  Figure  4  and  was 
first  demonstrated  by  Zoeppritz  in  1919.  According  to  Snell’s  Law,  for  the 
case  of  normal-incidence  compressional  (P-wave)  propagation  across  the  boun¬ 
dary  of  a  horizontally  oriented  system,  the  Zoeppritz  equation  for  the  ampli¬ 
tude  coefficients  of  the  incident  (Aj) ,  reflected  (A^),  and  transmitted  (A,)  waves 
becomes,  by  continuity  of  displacement,  simply 

A,-  =  A,  (3) 


and  for  continuity  of  stress  becomes 


E2/V2 

^1/^1 


Af 


(4) 


where  and  £2  ^^re  the  elastic  moduli  of  medium  1  and  2,  respectively,  and 
Vj  and  V2  are  the  velocities  in  each.  For  a  perfectly  elastic  medium  E  =  pv-, 
where  p  is  the  mass  density  and  v  the  elastic  P-wave  velocity,  the  quantity  pv 
is  called  the  acoustic  impedance  Z  of  the  medium  and  thus  represents  the 
influence  of  the  medium’s  characteristics  on  reflected  and  transmitted  waves. 
Solving  the  equations  of  stress  and  displacement  continuity  for  A^.  and  A,  in 
tenns  of  A-^  gives  * 


A 


C 


^1  ^2 
_! _ tA 

Zj  +Z2 


=  RA; 


(5) 
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and 


2Z 


1 


A  f  — —I  r  t  I 

'  Z1.Z2  ' 


=  TA: 


(6) 


The  coefficients  R  and  T  are  called  the  reflection  and  transmission  coef¬ 
ficients,  respectively,  and  represent  the  actual  percentage  of  the  wave’s 
reflected  and  transmitted  energy.  In  terms  of  the  reflection  coefficient,  the 
Zoeppritz  equations  can  be  written  as 


Z2 


\l-R) 


(7) 


providing  a  straightforward  method  for  determining  acoustic  impedance.  By 
knowing  the  first  Z  and  the  succeeding  Rs,  all  the  acoustic  impedances  can 
then  be  calculated.  In  this  case  the  first  layer  is  always  seawater,  which  has  a 
known  typical  impedance  value  of  1,550  lO^g/cm^s.  By  calculating  the 
remaining  Rs,  the  problem  is  solved.  However,  as  wiU  be  shown,  this  solution 
is  not  as  simple  as  it  seems. 


Classical  multilayer  reflectivity 

The  previous  discussion  presents  the  mechanics  of  seismic  wave  propaga¬ 
tion  across  a  single  boundary  between  two  media  without  accounting  for 
absorption  effects.  What  follows  is  the  case  of  a  multilayered  system  consist¬ 
ing  of  more  than  two  layers.  This  is  not  presented  as  a  derivation  of  the  math¬ 
ematics  involved:  rather,  it  is  provided  to  familiarize  the  reader  with  potential 
complexities  to  be  faced  in  any  attempt  to  develop  a  system  capable  of  eval¬ 
uating  multilayered  systems.  Excellent  detailed  presentations  of  the  mathemat¬ 
ics  of  propagation  in  layered  systems  are  found  in  the  work  of  Officer  (1958), 
Robinson  and  Treitel  (1980),  and  T.  M.  McGee  (1991). 

Discussion  will  be  restricted  to  compressional  (P-wave)  propagation  along 
the  vertical  axis  in  a  horizontally  layered  system.  For  convenience,  however, 
the  diagrams  are  drawn  with  time  along  the  horizontal  axis  so  that  the  rays 
appear  to  be  angled  through  the  layers.  For  this  discussion,  the  layers  are 
assumed  to  be  equal  travel-time  thicknesses.  Figure  5,  adapted  from  Robinson 
and  Treitel  (1980),  shows  the  ray  paths  of  the  multilayered  reflectivity  model, 
which  should  not  be  confused  with  the  commonly  seen  angled  ray-path  plots 
in  multilayered  systems  with  increasing  velocities. 

The  multilayered  system  presents  us  with  a  number  of  new  variables  not 
evident  in  the  simple  case  of  wave  propagation  across  a  single  boundary.  The 
effects  of  upward-traveling  waves  and  absorption  must  be  considered.  The 
relationship  between  wave  coefficients  is  defined  in  tenns  of  the  reflection  and 
transmission  coefficients.  Figure  6  presents  a  closeup  look  at  the  interface 
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between  layers  k  and  /:  +  1  at  time  /.  The  symbols  u  and  d  represent  upward- 
and  downward-traveling  waves,  respectively,  and  the  prime  symbol  signifies 
the  wave  is  at  the  bottom  of  the  layer.  Therefore,  the  coefficient  h’^/V  is  the 
upward-traveling  wave  in  the  bottom  of  layer  k,  comprising  the  part  due  to  the 
reflection  of  d j^i)  and  the  part  due  to  the  transmission  of  d as  follows: 

u[(i)  =  Rdl(i)*T^Uk^i(i)  (8) 

Likewise,  the  coefficient  is  given  as 

dk+l(i)  ^  Td\(i)+R' (9) 

These  equations  can  be  further  manipulated  to  yield  the  following  two 
simultaneous  equations  describing  the  situation  at  the  interface  at  time  /: 

Td^kd)  =  dk^j(i)^Ru^^j(i)  (10) 

and 

Tu\(i)  =  Rdi^^j(i)+u,^^j(i)  (H) 


Solution  of  the  simultaneous  equations  can  then  be  accomplished  using 
z-transform  theory  giving  the  expression  for  the  absorption-free  case  of  propa¬ 
gation  across  interface  i  in  the  matrix  form 


1 

'1 

1 

.D'tW- 

Rt  1 

V  J 

The  situation  becomes  further  complicated  when  the  effects  of  absorption 
are  introduced.  Absorption  is  a  frequency-dependent  attenuation  of  the  ampli¬ 
tude  of  a  waveform  as  it  propagates  through  a  given  sediment  layer.  In  order 
to  accurately  determine  the  reflection  coefficient —  and  thus  impedance —  at  a 
given  interface,  the  energy  loss  resulting  from  absorption  must  be  accounted 
for.  This  effect  can  be  modeled  in  the  above  approach  by  applying  the  appro¬ 
priate  absorption  operator.  A  proposed  solution  is  presented  by  T.  M.  McGee 
(1991)  and  will  not  be  repeated  here.  Specifics  of  the  AI  method,  offering  a 
practical  approach  to  the  problem,  will  be  presented  in  detail  later  in  this 
chapter. 
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Absorption  and  other  losses 

One  of  the  primary  energy  losses  encountered  during  acoustic  wave  propa¬ 
gation  through  differing  media  is  that  due  to  absorption.  Acoustic  reflections 
are  generated  from  impedance  mismatches  within  the  sediment  body.  The 
amount  of  returned  energy  depends  on  the  length  of  the  sound  path  and  the 
attenuation  of  sound  along  the  path.  Attenuation  along  the  sound  path  results 
from  a  number  of  mechanisms  as  described  by  Hamilton  (1980)  and  Sheriff 
(1975),  some  of  which  include: 

a.  Spherical  spreading  or  transmission  loss.  This  loss,  as  described  by 
Urick  (1983),  is  a  geometrical  effect  representing  the  regular  weakening 
of  a  sound  signal  as  it  spreads  outward  from  the  source.  As  will  be 
shown  later,  spreading  loss  varies  with  range  according  to  the  logarithm 
of  the  range. 

b.  Transmission  through  reflectors.  These  losses  include  multiple  reflec¬ 
tions,  reflection  and  refraction,  and  conversion  of  compressional  to 
shear  waves. 

c.  Reflector  roughness  and  curvature.  These  are  the  focusing  and  defo- 
cusing  effects  of  concave  and  convex  reflectors. 

d.  Scattering  due  to  inhomogeneities  within  the  sediment  body. 

e.  Intrinsic  absorption.  This  loss  involves  a  process  of  conversion  of 
acoustic  energy  into  heat  and  thereby  represents  a  true  loss  of  acoustic 
energy  to  the  medium  in  which  propagation  is  taking  place  (see 
Equation  14). 

A  discussion  concerning  the  assessment  of  absorption  will  be  addressed 
later  in  the  sections  “Absorption  Calibration”  and  “Spectrum  Analysis.”  The 
total  of  all  losses  is  called  the  effective  attenuation,  and  because  of  the  ran¬ 
domness  of  items  c  and  d  above  (and  many  other  possible  factors  not  men¬ 
tioned),  it  is  nearly  impossible  to  completely  accoimt  for  in  the  real  world  of 
acoustic  profiling.  However,  the  major  sources  of  attenuation  (i.e.,  spreading, 
reflection,  and  absorption)  have  been  researched  extensively,  providing  reason¬ 
able  approximations  of  the  actual  losses  occurring.  For  the  case  of  absorption 
in  marine  sediments,  there  has  been  considerable  debate  concerning  the  most 
appropriate  attenuation  model.  Stoll  (1974)  following  Biot  (1956a  and  1956b) 
favore  a  theoretical  model  where  at  low  frequencies  (generally  less  than 
10  kHz)  absorption  has  an/^  dependence  and  at  the  higher  frequencies  an  / 
dependence.  Hamilton  (1972a  and  1972b)  and  others  (refer  to  references  in 
Hamilton  papers),  on  the  other  hand,  present  convincing  experimental  evidence 
as  to  absorption’s  relationship  to  the  first  power  of  frequency.  Hamilton 
(1972a)  presents  the  following  important  observations: 

a.  Absorption  is  dependent  on  the  first  power  of  frequency. 
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b.  Velocity  dispersion  is  insignificant  when  compared  to  a. 

c.  Intergrain  friction  appears  to  be,  by  far,  the  dominant  cause  of  wave- 
energy  dampening  in  marine  sediments. 

Specifically,  absorption  varies  as  a  function  of  frequency  according  to  the 
empirical  equation 

fy  =  w"  (13) 


where 


a  =  absorption,  db/m 
f  =  frequency,  kHz 
k  =  attenuation  coefficient,  dbImIkHz 
n  =  exponent  of  frequency 

The  constant  n  has  been  experimentally  determined  to  be  essentially  unity 
for  the  frequencies  of  interest  figure  7)  leaving  k  in  Equation  13  as  the  only 
variable.  This  constant  varies  with  sediment  type  and  is  related  to  porosity 
and  mean  grain  size  as  shown  in  Figure  8.  Because  of  the  extensive  experi¬ 
mental  data  utilized,  Hamilton’s  linear  relationship  was  chosen  as  the  model 
for  the  AI  method.  A  modification  of  this  model  as  described  by  Caulfield 
and  Yim  (1983)  and  Caulfield,  Caulfield,  and  Yim  (1985)  is  utilized  in  the  AI 
method  to  estimate  the  engineering  properties  of  marine  sediments.  A  reason¬ 
able  measure  of  absorption,  in  keeping  with  Equation  13,  is  provided  assuming 
an  exponential  correction  as  a  fimction  of  frequency  by 


a  =  10  logioc 


(14) 


where 

p  =  density  of  layer,  gm/cc 
c  =  soimd  velocity  of  layer 
k  =  attenuation  coefficient  (similar  to  Hamilton’s) 

X  =  precision  absorption  correction  factor 

The  coefficient  k  is  either  experimentally  derived  or  estimated  from  Hamil¬ 
ton’s  regression  equation  (refer  to  Figure  8).  A  correction  factor  X  is  included 
to  compensate  for  localized  variations  in  the  absorption  properties  of  sediments 
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in  a  given  geologic  setting.  This  value  normally  remains  unity,  is  applied  in  a 
linear  manner,  and  is  only  altered  when  detailed  core  data  are  available.  The 
value  of  X  is  increased  or  decreased  so  that  the  acoustically  derived  impedance 
estimates  match  the  laboratory  measured  core  properties.  It  is  important  to 
note  that  research  regarding  the  acoustic  absorption  in  marine  sediments  has 
been  widely  studied.  A  compilation  of  the  various  absorption  theories  and 
models  is  provided  by  Hampton  (1974).  The  model  utilized  in  the  A1  method 
is  by  no  means  totally  perfect  in  terms  of  applicability  to  all  possible  marine 
environments.  Rather,  it  is  designed  to  provide  a  reasonable  estimate  towards 
the  prediction  of  actual  sediment  properties  and,  upon  critical  examination  of 
actual  in  situ  conditions,  to  be  refined  to  precisely  model  a  particular  sediment 
environment. 


Relationship  of  Acoustic  Impedance  to 
Geotechnical  Properties 

Because  acoustic  impedance  basically  represents  the  influence  of  a  medi¬ 
um’s  characteristics  on  reflected  and  transmitted  waves,  many  geotechnical 
properties  such  as  porosity,  density,  mean  grain  size,  bulk  modulus,  etc.,  exhi¬ 
bit  excellent  correlation  with  impedance.  During  the  last  two  decades,  the 
ability  to  predict  geotechnical  properties  from  nonnal  reflectivity  through 
impedance  calculations  has  become  very  much  established.  However,  seismic 
signatures,  and  therefore  acoustic  impedances,  are  not  considered  unique.  Sev¬ 
eral  combinations  of  geologic  conditions  could  conceivably  yield  similar  signal 
characteristics.  Therefore,  a  critical  stage  in  the  process  is  development  of 
geoacoustic  relationships  which  are  used  to  model  a  specific  geologic  environ¬ 
ment.  General  relationships  for  acoustic  impedance  versus  density  are  given 
by  Hamilton  and  Bachman  (1982)  in  Figure  9  and  for  soil  type  by  Caulfield 
and  Yim  (1983)  in  Table  1.  These  relationships  (or  estimates)  are  based  on 
worldwide  averages  of  impedance  versus  sediment  properties  from  primarily 
deep  ocean  environments  and  do  not  necessarily  constimte  the  precise  charac¬ 
terization  of  all  geologic  situations. 

The  AI  method  utilizes  a  regional  calibration  approach.  For  a  particular 
study  area  (for  example,  Panama  City  Harbor  or  Delaware  Bay),  computed 
acoustic  impedance  values  are  related  directly  to  the  soil  parameters  measured 
from  at  least  ten  cores  representative  of  the  area.  A  recently  completed  study 
in  the  predominantly  sand  environment  off  the  coast  of  Delaware  by  McGee 
(in  publication)  presents  an  excellent  example  of  the  regional  approach  to  geo¬ 
acoustic  modeling.  Figure  10  presents  the  relationship  between  computed 
impedance  and  the  soil  properties  of  density  and  mean  grain  size  for  the  Dela¬ 
ware  Coast  project.  It  is  important  to  note  that  for  the  Delaware  model,  the 
predictions  are  only  applicable  for  the  sediment  ranges  presented.  Figure  1 1 
shows  an  example  of  density  versus  mean  grain  size  derived  from  laboratory 
measurements  of  core  samples  collected  by  the  U.S.  Anny  Engineer  Water¬ 
ways  Experiment  Station  (WES)  from  Panama  City  Harbor  and  Gulfport  Ship 
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Table  1 

Soil  Classification  versus  Acoustic  impedance 

Description 

Acoustic  Impedance 

iqpf.g 

cm^  s 

Water 

1450  -  1550 

Silty  Clay 

2016  -  2460 

Clayey  Silt 

2460  -  2864 

Silty  Sand 

2864  -  3052 

Very  Fine  Sand 

3052  -  3219 

Fine  Sand 

3219-3281 

Medium  Sand 

’  3281  -  3492 

Coarse  Sand 

3492  -  3647 

Gravelly  Sand 

3647  -  3880 

Sandy  Gravel 

3880  -  3927 

Channel.  Core  samples  from  Savannah  Ship  Channel  were  obtained  from  the 
South  Atlantic  Division  Laboratory. 

These  examples  by  no  means  constitute  the  entire  library  of  impedance 
versus  soil  properties  accumulated  during  development  of  the  AI  method.  As 
a  matter  of  fact,  geoacoustic  relationships  are  under  continual  refinement  as 
new  surveys  are  performed  in  geologically  different  areas.  Geoacoustic  data 
have  been  collected  in  over  17  areas  to  date.  These  include  work  performed  in 
Pacific  coast  harbors,  the  Gulf  of  Mexico,  Atlantic  coastal  harbors,  and  some 
inner  rivers  and  reservoirs.  More  specifics  on  these  areas,  including  a  com¬ 
plete  listing,  will  be  presented  in  Chapter  4,  “Case  Studies.”  Research  utilizing 
this  extensive  database  is  ongoing  at  WES  in  the  development  of  comprehen¬ 
sive  geoacoustic  relationships  of  sediments  representative  to  the  riverine,  har¬ 
bor,  and  near-coastal  environments  typical  of  Corps  projects.  A  statistical 
analysis  of  acoustic  data  collected  from  several  of  the  above  sites  is  presented 
at  the  conclusion  of  Chapter  4,  “Case  Studies.” 


Engineering  Approach 

Utilizing  proper  calibration  procedures  with  data  of  high  signal-to-noise 
ratio,  seismic  reflection  data  can  be  processed  to  accurately  estimate  the  den¬ 
sity  and  soil  type  of  bottom  and  subbottom  sediments.  Site-specific  calibra¬ 
tions  are  perfonned  on  every  Job  by  correlating  acoustic  impedance  calculated 
from  seismic  reflection  data  at  a  core  location  with  in  situ  infonnation  (den¬ 
sity,  mean  grain  size,  etc.)  at  that  location.  Experience  to  date  has  shown  that 
calibrations  made  at  a  few  locations  within  a  geologic  region  produces  the 
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shallow  seismic  parameters  necessary  to  adequately  calibrate  and  describe  the 
entire  region  as  presented  in  R.  G.  McGee  (1991);  Ballard  and  McGee  (1991); 
Ballard,  McGee,  and  Harmon  (1991);  Ballard,  McGee,  and  Whalin  (1992); 
Ballard  et  al.  (1992);  Ballard  et  al.  (1993);  Sjostrom,  Ballard,  and  McGee 
(1991,  1992);  Sjostrom,  McGee,  and  Ballard  (1992). 


Source  calibration  model  and  reflection  coefficient 

Sonar  equations.  The  geoacoustic  parameter  calibration  procedure  begins  by 
determining  the  total  acoustic  energy  incident  at  the  bottom  surface.  This 
basically  involves  determining  the  precise  reflection  coefficient  for  the  first 
reflector  (bottom  surface)  and  its  associated  acoustic  bottom  loss  for  a  given 
sediment.  Since  the  sound  velocity  of  water  and  its  density  can  be  readily 
measured,  the  absolute  impedance  of  water  can  be  calculated.  Knowledge  of 
the  reflection  coefficient  for  the  water  bottom  interface,  which  is  completely 
independent  of  frequency,  allows  direct  computation  of  the  absolute  impedance 
of  the  first  layer  of  the  bottom.  The  total  energy  produced  by  the  source,  or 
source  wavelet,  must  be  known  absolutely.  This  is  accomplished  through  use 
of  a  calibration  hydrophone  allowing  determination  of  source  level  (SL)  and 
transmission  losses  associated  with  underwater  acoustic  wave  propagation 
through  the  sonar  equations.  The  sonar  equations,  discussed  thoroughly  by 
Urick  (1983),  describe  quantitative  effects  on  sonar  equipment  created  by 
many  phenomena  peculiar  to  underwater  sound  production.  These  equations 
are  both  design  and  prediction  tools  for  underwater  sound  applications  and 
relate  effects  of  the  medium,  target,  and  equipment.  The  general  sonar  equa¬ 
tion  is  given  as  follows: 

Sf^  =  SL-N^-N^y^^N^^DI*BL  (15) 

where 

Sf^  =  bottom  reflection  energy  at  receiver,  db 

SL  =  total  energy  of  source,  db 

=20  X  LogjQ  (range,  meters),  db  (transmission  loss  due  to  spherical 
spreading  along  the  path  of  propagation) 

^hyd  ~  receiver  sensitivity,  db 
=  amplifier  gain,  db 

DI  =  directivity  index  of  receiving  array,  db  (function  of  transducer  beam 
pattern) 

BL  =  bottom  loss,  db  =  20 
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R  =  reflection  coefficient 


The  following  sections  of  this  report  will  discuss  the  many  variations  and 
utility  of  the  sonar  equation  in  perfonning  precise  field  calibrations  of  equip¬ 
ment  and  surface  reflections.  Figure  12  is  a  detailed  depiction  of  tlie  physical 
elements  in  a  normal  calibration  and  bottom  reflection  sonar  equation  solution 
case.  The  value  includes  all  preamplifiers  and  amplifiers  and  is  obtained 
from  the  electrical  calibration  of  the  receiving  equipment  Figure  13  shows  a 
typical  amplifier  gain  calibration  curve  necessary  for  determination  of  absolute 
energy.  The  calibration  hydrophone  sensitivity  available  from  man¬ 

ufacturers  of  the  hydrophone  and  should  be  traced  to  the  ANSI  Standard.  The 
receiving  array  sensitivity  may  also  be  available  from  the  manufacturer 

or  can  be  easily  calibrated  in  the  field  using  the  calibration  hydrophone  and  an 
alternate  form  of  the  sonar  equation.  This  procedure  will  be  discussed  in 
detail  later  in  this  report. 

A  suite  of  signal  analysis  programs  has  been  developed  by  WES  in  the 
course  of  this  study  specifically  to  perfonn  all  aspects  of  seismic  calibrations. 
Additionally,  these  programs  are  used  to  evaluate  all  parameters  of  the  sonar 
equations  and  to  analyze  measured  acoustic  signals  representative  of  bottom 
conditions. 

Directivity  index.  The  directivity  index  (DI)  is  a  function  of  the  beam  pat¬ 
tern  of  the  transducer  array  and  is  an  indication  of  the  amount  of  the  total 
signal  the  hydrophone  is  pennitted  by  its  sensitivity  pattern  as  shown  in  Fig¬ 
ure  14.  The  higher  the  DI,  the  more  discriminating  the  hydrophone  is  against 
signals  arriving  from  directions  other  than  along  the  desired  acoustic  axis. 
Figure  14  presents  the  directional  pattern  of  the  MASSA  Model  TR-75A  trans¬ 
ducer,  typical  of  transducers  utilized  in  the  “pinger”  class  of  profiling  systems. 
Because  WES  typically  offsets  its  receivers  horizontally  from  the  transmitter, 
as  explained  in  Appendix  A,  the  DI  becomes  a  significant  parameter  due  to 
reflection  angles  along  the  path  of  propagation.  Figure  15  presents  (a)  equip¬ 
ment  geometry  for  the  Delaware  coast  survey  and  its  effect  on  directivity,  and 
(b)  the  DI  correction  versus  water  depth  for  application  in  the  sonar  equation. 

Source  level  (SL)  calibration.  The  first  step  in  the  calibration  process  is  to 
determine  the  absolute  source  level.  Some  sonars  have  this  data  available  from 
the  manufacturers.  Unfortunately,  many  seismic  systems  do  not  have  this 
information  readily  available,  and  even  if  they  did,  field  operating  conditions 
vary  to  such  an  extent  that  the  published  levels  are  not  sufficient  for  precise 
reflection  computations. 

The  direct  wave  calibration  of  the  sonar  source  level  is  accomplished  by 
writing  the  sonar  equation  for  the  measurement  of  the  direct  wave  via  a  cali¬ 
bration  hydrophone  as  follows: 
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where 


=  direct  wave  signal  level,  db 

^wdir  ~  transmission  loss  between  source  and  cal  phone,  db 

All  the  terms  in  Equation  16,  except  SL,  are  either  absolutely  known  or 
directly  measured.  Therefore,  solving  for  SL,  the  absolute  source  level  for  any 
particular  seismic  system  can  be  determined. 

The  calibration  hydrophone  is  introduced  into  the  sound  field  along  the 
acoustic  axis  of  the  source.  The  hydrophone  must  be  no  closer  than  about  5  m 
to  the  source  to  avoid  near-field  wave  form  problems.  Normally  calibration 
data  arc  taken  at  various  depths  to  vary  the  transmission  loss  and  to  confirm 
that  the  hydrophone  is  veiy  near  the  acoustic  axis.  The  calibration  signal 
should  obey  the  inverse  square  law;  i.e.,  20  logjQ  (ej  -  62)  =  20  logjQ  (r2  -  rj) 
where  Cj  is  the  voltage  generated  by  the  hydrophone  at  range  and  62  is  the 
voltage  generated  by  the  hydrophone  at  range  r2. 

Figure  16  presents  a  typical  seismic  system  calibration  data  plot  recorded 
during  a  recent  field  exercise.  This  single  data  record  contains  all  the  field 
data  required  to  completely  calibrate  all  aspects  of  equipment  operations  and 
provides  calibration  data  for  the  surface  sediment  impedance.  The  source  level 
calibration  is  performed  utilizing  the  data  between  file  number  0002  and  0004 
where  variations  in  amplifier  gain  and  hydrophone  range  are  occurring.  An 
example  calculation  is  provided  in  Figure  17  on  data  collected  in  Delaware 
Bay  just  off  Cape  May  Inlet,  NJ.  The  figure  presents  results  spatially  averaged 
over  40  con^cutive  soundings  with  the  table  showing  the  statistical  variations 
typical  of  acoustic  profiling.  The  statistical  properties  of  the  SL  (source  level) 
and  BL  (bottom  loss)  are  given  in  decibels  for  the  subfile  pci  10012  as  well  as 
for  the  sequence  of  files  between  subfiles  0012  and  0015. 

Receiving  hydrophone  sensitivity  calibration.  As  with  the  source  level,  the 
array  sensitivity  of  receiving  hydrophones  must  be  absolutely  known. 

Field  calibration  is  performed  by  comparing  the  signal  levels  of  the  receiving 
array  with  the  calibration  hydrophone  over  the  same  bottom  condition.  The 
calibration  hydrophone  is  located  in  the  immediate  vicinity  of  the  receiving 
array  at  the  same  depth  elevation  as  presented  in  Figure  18.  The  sonar  equa¬ 
tion  is  designed  to  solve  for  Ni,ydr  follows: 

t^hydr  =  %dc^SprSRc-NAr-^Ac 


where  S^^  and  S^^,  are  the  received  signals  and  amplifier  gains  for 
the  receiving  array  and  calibration  hydrophone,  respectively.  Figure  19  is  an 
example  calculation  from  the  same  calibration  trial  discussed  for  SL 
determination. 
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Determination  of  bottom  loss  and  surface  reflection  coefficient.  Once  all 
equipment  parameters  are  known,  the  transmission  loss  defined,  and  all  deploy¬ 
ment  parameters  (such  as  Dl)  assessed,  the  bottom  surface  characteristics  can 
be  ev^uated  through  the  sonar  equation  by  rearranging  Equation  15  to  solve 
for  BL  as  follows: 

BL  =  SR^N^yj-SL^N^-N^-Dl  (18) 

Since  all  terms  on  the  right  side  of  the  equation  are  now  known,  BL,  and 
therefore  the  surface  reflection  coefficient  (BL  =  20  logjfJR)  and  acoustic 
impedance  (Equation  7),  can  be  readily  determined  as  shown  in  Figure  17.  If 
the  desired  result  is  an  assessment  of  the  bottom  surface  characteristics,  the 
acoustic  solution  is  complete.  AH  that  remains  is  the  correlation  of  the  acous¬ 
tic  parameters  with  physical  sediment  properties  through  geoacoustic  relation¬ 
ships  (Table  1,  Figures  9  and  10).  Figure  20  shows  a  typical  surface  sediment 
characterization  in  terms  of  mean  grain  size  and  density. 

Another  straightforward  method  for  determining  the  acoustic  bottom  loss  is 
the  utilization  of  multiple  reflections  when  present.  The  difference  in  energy 
between  the  bottom  echo  wavelet  and  the  multiple  wavelet  over  the  same 
bottom  is  simply  the  sum  of  transmission  loss  and  bottom  loss.  Figure  21  pre¬ 
sents  this  relationship  in  terms  of  the  sonar  equation  with  an  example  calcula¬ 
tion.  It  is  important  to  note  that  this  method  is  only  good  when  high-quality 
multiple  data  are  available.  This  is  a  function  of  the  frequency  of  the  sound 
source  and  the  relative  competence  of  bottom  materials.  Also,  survey  condi¬ 
tions,  such  as  a  high  sea  state,  or  equipment  deployment  can  affect  the  quality 
of  multiple  data. 

As  stated  earlier,  with  the  first  Z  and  the  succeeding  Rs  known,  aU  the 
acoustic  impedances  can  be  calculated.  All  that  remains  is  to  apply  the  appro¬ 
priate  absorption  model  to  the  acoustic  data. 


Subbottom  reflection  coefficients  (absorption  calibration) 

An  overview  of  the  model  utilized  as  an  initial  estimate  of  the  absorption 
characteristics  of  marine  sediments  has  been  provided  by  Equations  13  and  14. 
Caulfield  and  Yim  (1983)  describe  this  in  some  detail.  There  are  a  number  of 
theories  on  the  absorptive  properties  of  different  marine  sediments,  any  of 
which  might  be  better  suited  to  a  specific  site  than  the  Caulfield  model. 
Therefore,  the  procedure  allows  the  use  of  different  absorption  algorithms, 
basically  through  modification  of  the  velocity  versus  density  relationships 
utilized  in  the  AC50  program  (Caulfield  1992a).  However,  this  should  only  be 
done  when  corroborating  core  data  are  available  to  verify  the  model  prediction. 

The  software  procedure  was  developed  (Caulfield  1992a)  to  provide  both 
“precision”  analysis  of  unique  seismic  traces  and  “production-type”  processing 
of  continuous  seismic  profile  lines.  The  AC50  program,  developed  through 
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this  work  unit,  is  described  in  detail  in  the  above  reference  and  basically  fol¬ 
lows  the  methodology  described.  Once  SL  calibrations  and  initial  BL  charac¬ 
teristics  are  determined,  providing  initial  calibration  coefficients,  the  AC50 
program  utilizes  the  classical  multilayer  reflectivity  algorithm  with  appropriate 
absorption  operator  to  determine  reflection  coefficients  and  associated  acoustic 
impedances  for  each  succeeding  reflection  horizon. 

Specifically,  each  data  trace  is  converted  to  an  energy  matrix  yielding 
energy  versus  travel  time.  Since  the  program  now  is  designed  to  handle  only 
the  band-limited  center  fi'equency  case,  this  energy  matrix  consists  of  only  one 
vector  of  energy  values  dt^)  where  (o^  is  the  center  frequency  and  dt^ 
is  the  data  window  subsection.  Each  matrix  element  is  then  corrected  for 
absorption.  An  initial  assumption  is  made  concerning  the  average  density  and 
sediment  velocity  based  on  the  duration  of  the  bottom  signal.  A  long-duration 
signal  usually  implies  high  clay  content,  a  short  signal  a  predominance  of  sand, 
Energy  matrix  elements  corrected  for  absorption  are  in  the  form 


E'e 


P®«i 

kc 


(19) 


where  the  exponent  is  the  same  as  in  Equation  13.  Now  that  the  initial 
absorption  correction  has  been  applied  for  the  frequency  utilized,  the  total 
energy  per  layer  can  be  calculated  as 


M 

m=l 


(20) 


where 

Ef  =  total  energy  per  layer 

m  =  sample  number 

M  =  total  number  of  samples 

^corr  =  ^abs  +  Spherical  spreading  correction 

The  reflection  coefficient  of  the  layer  is  the  ratio  of  the  energy  from  that 
layer  divided  by  the  total  energy  incident  on  the  bottom  minus  the  energy 
reflected  from  the  previous  layers.  After  the  reflection  coefficients  have  been 
computed,  the  sediment  impedance  for  any  layer  may  be  calculated  by 
Equation  7. 

If  core  data  are  available,  this  process  is  repeated,  changing  the  absorption 
correction  factor  until  the  layer  estimates  match  those  of  the  core  data.  Sedi¬ 
ment  characterization  is  then  accomplished  using  the  appropriate  geoacoustic 
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model  and  should  apply  for  the  site  under  investigation.  Examples  of  outputs 
generated  during  this  process  are  given  in  Figures  22-25.  A  detailed  study  of 
each  parameter  in  the  process,  including  sensitivity  analysis  of  each  of  the 
variables,  is  provided  by  Caulfield  (1990).  Ground  truth  examples  are  pro¬ 
vided  in  Chapter  4,  “Case  Studies.” 


Spectrum  analysis 

Knowledge  of  the  spectral  content  of  acoustic  signals  (frequency  and  total 
energy)  is  vital  to  the  accurate  application  of  the  process.  Frequency  content 
of  a  seismic  trace  is  highly  useful  because  the  frequency  content  governs  both 
temporal  and  spatial  resolution.  The  spectral  response  of  any  portion  of  the 
data  also  indicates  the  absorptive  nature  of  overlying  sediments,  since  different 
sediments  have  characteristic  absorptive  tendencies.  Also,  as  stated  in  previous 
sections,  not  all  of  the  energy  loss  exhibited  in  the  data  is  due  to  absoiption. 
The  energy  loss  beneath  a  hard  layer,  for  instance,  is  due  both  to  classical 
absorption  and  to  the  fact  that  a  large  percentage  of  the  energy  is  reflected 
upwards.  Since  higher  frequencies  resolve  thin  beds  better,  it  may  be  that  a 
thin,  high  impedance  bed  will  selectively  reflect  higher  frequency  seismic 
energy.  The  spectral  analysis  will  provide  insight  into  the  absorptive  prop¬ 
erties  of  unique  sediments  by  measuring  the  distribution  of  energy  content 
across  the  spectrum,  leading  to  the  classification  of  the  sediments. 

Spectral  analysis  is  also  quite  useful  in  differentiating  random  noise  from 
true  seismic  reflections.  By  studying  the  effects  of  different  filters  and  the 
signatures  of  various  seismic  sources,  data  acquisition  can  be  optimized. 

The  Digital  Spectral  Analysis  (DSAIO)  shallow  seismic  software  (Caulfield 
1991b)  was  developed  to  provide  seismic  amplitude  and  power  spectra  versus 
time  analysis  using  Fourier  transforms  of  the  DFAS  digital  data.  Figures  26  - 
28  show  typical  time-frequency  plots  produced  by  the  DSAIO  program  for 
3.5-kHz  (pinger)  and  boomer  reflection  data,  and  chirp  source  wavelet  data. 

On  the  far  left  of  the  figures  is  the  amplitude  record  for  a  single  subfile  of  dig¬ 
ital  data,  which  is  a  collection  of  40  consecutive  soundings.  A  single  ampli¬ 
tude  trace  (usually  the  last  sounding  in  the  subfile)  is  plotted  with  signal 
amplitude  in  volts  plotted  along  the  horizontal  axis.  The  vertical  axis  for  all 
plots  is  time  and  represents  depth.  The  remainder  of  the  display  is  the  spectral 
analysis  plot.  The  horizontal  axis  represents  frequency,  increasing  to  the  right. 
Each  colored  box  represents  the  energy  at  a  particular  frequency  in  time.  A 
legend  relating  color  to  energy  levels  is  provided  on  Figure  26.  The  darker  the 
color,  the  higher  the  energy  content  of  the  frequency  component.  To  the  right 
of  the  spectral  plot  are  the  parameters  used  in  the  display.  An  ASCII  file  (not 
shown)  of  all  the  time-frequency  energy  values  provides  detailed  analysis. 

This  display  has  a  dynamic  range  of  over  60  db,  and  consequently  the  “empty” 
water  column  may  have  frequency  values  corresponding  to  random  noise. 

By  comparing  the  energies  at  specific  frequencies  between  reflection  hori¬ 
zons  and  compensating  for  reflection,  energy  loss  due  to  absorption  can  be 
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determined.  Also,  by  evaluating  attenuation  across  the  spectrum  between 
horizons,  layer  characteristics  can  be  estimated.  The  DSAIO  program  is  an 
important  tool  required  to  develop  and  apply  the  appropriate  absorption  operat¬ 
ors  for  the  subbottom  (absorption)  calibration  in  the  AC50  program. 


Layer  Identification  Methods 

The  procedure  described  so  far  computes  impedances  as  a  function  of  the 
reflection  coefficient  at  layer  boundaries.  In  the  purely  classical  model  of 
natural  marine  sediments,  the  lithostratigraphy  is  usually  modeled  as  multiple 
layers  of  increasing  impedances,  and  for  this  case  the  model  presented  here 
works  beautifully.  However,  as  stated  earlier,  several  combinations  of  geologic 
conditiCMis  could  exhibit  very  similar  signal  characteristics,  resulting  in  incor¬ 
rect  sediment  characterizations  based  solely  on  acoustic  impedance.  Therefore, 
a  number  of  different  tools  have  been  developed  to  assist  the  engineer  in  cor¬ 
rectly  assessing  the  geoacoustic  model. 

Correlation  processing.  Correlation  analysis  aids  in  detennining  the  gen¬ 
eral  statistical  properties  of  sonar  signals,  allowing  optimization  of  the 
detection  process  by  improving  both  the  signal-to-noise  (S/N)  and  the  time 
resolution.  Also,  the  correlation  of  good  S/N  broad-band  data,  such  as 
“chiiped”  signals,  provides  information  for  detennining  the  polarity  of  the 
reflection  coefficient.  The  Digital  Shallow  Seismic  Processing  and  Correlation 
System  (DPCIO),  developed  to  meet  these  objectives,  is  generally  discussed 
below  and  is  described  in  detail  by  Caulfield  (1991c). 

The  correlation  process  employed  is  a  matched  filter,  equivalent  to  cross- 
correlation  of  the  received  reflection  data  with  the  source  wavelet.  In  other 
words,  if  two  wavelets  are  identical,  the  filter  output  is  simply  the  auto-cor¬ 
relation  function  of  the  source  wavelet,  and  the  filter  is  said  to  be  “matched”  to 
the  source  wavelet.  The  cross-correlation  function  of  a  time  series  X(t)  with 
the  source  wave  form  Y(t)  is: 


X„(0*X„(r) 


(21) 


where 

X^(t)  =  amplitude  of  the  sample 

T  =  lag  or  offset  in  Ar  (normally  the  sampling  interval) 

In  this  version  of  the  correlation  process,  the  reference  signal  is  provided  as 
an  rt-sample  length  source  signal  wavelet,  typically  measured  directly  with  the 
calibration  hydrophone  as  described  in  the  source  level  calibration  section. 

This  signal  is  then  correlated  against  the  reflection  data  and  output  is  digitally 
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stored  in  the  DFAS  output  fonnat.  To  ensure  consistency  and  to  allow  AC50 
processing  of  the  correlated  data,  the  program  ensures  that  the  energy  content 
of  the  raw  signal  trace  and  the  correlated  signal  trace  are  equal.  The  data  are 
also  nomialized  to  maintain  the  relative  amplitudes  of  the  reflection  signals. 
Figures  29-31  show  the  improvements  in  S/N  and  resolution  possible  through 
correlation  analysis  of  noisy  pinger,  boomer,  and  chirped  signals.  Detailed 
examples  of  the  benefits  of  correlation  processing  on  a  wide  variety  of  seismic 
data  are  illustrated  in  Caulfield’s  manual  (1991c).  Figure  32  shows  the 
improvement  in  S/N  provided  by  correlation  processing  on  boomer  data  col¬ 
lected  during  a  previous  survey.  Since  this  correlation  process  is  utilized  in 
many  of  the  layer  identification  methods  described  below,  specific  examples  of 
the  correlation  applications  will  be  presented  under  tliat  topic. 

Synthetic  forward  modeling.  An  important  tool  utilized  in  layer  identifica¬ 
tion,  project  planning  and  feasibility  studies,  and  optimization  of  data  acquisi¬ 
tion  and  analysis  methodologies  is  the  process  of  manufacturing  synthetic 
seismograms  through  geoacoustic  forward  modeling.  The  modeling  techniques 
developed  for  this  method  involve  creating  one-dimensional  synthetic  seismo¬ 
grams  by  simply  convolving  a  known  wave  fonn,  such  as  the  source  wavelet 
measured  during  calibration,  with  a  reflectivity  Ihnction  derived  from  actual 
geologic  data.  The  simplest  and  therefore  most  straightforward  approach  is  the 
reflectivity  model  utilizing  primary  waves  only,  with  no  absorption  corrections, 
and  assuming  normal  incidence.  For  very  shallow  (upper  20  ft  of  subbottom), 
relatively  thin  sediment  units,  this  assumption  is  usually  sufficient.  The  mod¬ 
eling  technique  described  below  includes  correlation  processing  and  decon¬ 
volution  processing  and  is  discussed  as  an  example  calculation. 

Two  geologic  situations  very  common  to  the  near-shore  enviromnent  were 
considered:  (1)  sand  over  soft  mud  over  hard  clay,  and  (2)  sand  over  com¬ 
pacted  clay  (Figure  33).  The  properties  of  each  sediment  unit  were  modeled 
and  are  shown  with  each  figure  as  computed  impedances  and  reflection  coef¬ 
ficients.  Source  wavelets  for  a  boomer  signal,  a  shaded  1-  to  6-kHz  chirp,  and 
a  7-kHz  pinger  are  convolved  with  the  reflection  model  to  generate  simulated 
subbottom  reflection  data  as  shown  in  Figures  34  and  35.  The  chirp  data  must 
be  correlated  with  the  source  wavelet  before  it  provides  any  useful  infor¬ 
mation.  Whereas  the  seismograms  provide  reasonable  pictures  of  the  mag¬ 
nitudes  of  the  reflections  expected  at  each  horizon,  the  pinger  and  boomer  data 
show  vertical  resolutions  on  the  order  of  only  1 .5  ft  and  in  practical 
applications  make  it  impossible  to  identify  the  sign  (±)  of  the  reflector.  How¬ 
ever,  the  processes  of  correlation  and  deconvolution  can  assist  in  making  these 
detenninations. 

The  correlation  technique  is  used  to  perfomi  the  matched  filter  process 
described  previously  and,  for  the  boomer  and  chirp  data,  improves  botli  llie 
S/N  and  resolution.  An  example  of  the  correlated  boomer  signal  for  reflection 
sequences  is  provided  in  Figure  36  and  shows  a  factor-of-3  improvement  in 
resolution  and  an  indication  of  the  sign  of  the  reflection  coefficient. 


Chapter  2  Technical  Development 


Correlation  of  the  shaded  chirp  signal  for  the  reflection  sequences  provides 
improvement  in  vertical  resolution  by  more  than  a  factor  of  4,  as  shown  in 
Figure  36.  Events  within  0.5  ft  can  be  resolved.  Also,  some  indication  of 
reflection  sign  is  shown. 

The  process  of  deconvolution  involves  generating  a  filter  that  transforms 
the  calibrated  source  wavelet  into  a  spike.  Figure  37  illustrates  the  results  of 
deconvolution  as  applied  to  boomer  data  for  both  reflection  sequences.  It  can 
be  seen  that  the  resolution  of  peaks  can  be  easily  resolved  to  within  0.3  fL 
More  importantly,  the  sign  of  the  reflection  coefficient  is  readily  determined, 
allowing  for  easy  identification  of  sands  over  soft  clays.  The  same  spiking 
deconvolution  can  be  applied  to  correlated  chiip  data  for  further  improving 
resolution  and  sign  identification. 

The  modeling  process  is  also  utilized  to  assist  in  the  interpretation  of  sedi¬ 
ment  characteristics  by  comparing  modeled  results  with  actual  subbottom  data. 
An  example  of  this  approach  is  presented  in  Figure  38  for  boomer  data  col¬ 
lected  in  an  area  where  fine  sediments  were  bounded  above  and  below  by 
competent  sands  and  gravels.  Actual  boomer  data,  both  raw  and  correlated, 
are  compared  to  the  synthetic  seismogram  generated  from  the  reflection 
sequence  developed  from  core  log  data  obtained  from  the  surveyed  area.  The 
data  show  a  thin  layer  of  high-impedance  sand  overlying  clayey  sands  and  silty 
clays  overlying  poorly  graded  gravels.  Note  the  polarity  of  the  reflection 
coefficients  and  how  the  correlated  boomer  trace  provides  some  indication  of 
the  direction  of  impedance  change.  In  this  case,  the  modeling  process  pro¬ 
vided  furtlter  verification  of  the  geoacoustic  analysis. 

Polarity  of  reflection  co^icient.  Broad-banded  data  having  good  S/N,  such 
as  chirper  or  even  boomer  data,  can  be  analyzed  using  correlation  processing 
to  provide  an  indication  of  the  nature  of  the  impedance  change  (higher  or 
lower)  at  reflection  horizons.  For  example,  a  soft  layer  underneath  a  hard 
layer  will  produce  a  minus  impedance  change,  or  negative  reflection  coeffi¬ 
cient  This  negative  reflection  coefficient  results  from  the  phase  change  of  the 
reflected  signal  occurring  when  the  wave  reflects  off  a  softer  layer.  The 
match-filter  correlation  of  the  source  wavelet  with  the  reflected  wave  from  a 
deeper  soft  layer  results  in  a  “flipped”  or  negative  reflection  signal  as  shown 
by  Figure  38.  This  technique  works  only  on  data  with  good  S/N  and  is  very 
difficult,  if  not  impossible,  to  perform  on  band-limited  data  such  as  those  pro¬ 
duced  by  pingers.  A  beta  version  of  an  updated  AC50  program  incorporates 
estimation  of  the  reflection  coefficient  polarity  with  very  encouraging  results. 

At  the  time  of  this  writing,  however,  one  must  detennine  the  sign  (+  or  -) 
outside  the  AC50  program  through  either  correlation  processing  or  determina¬ 
tion  of  the  signal  phase.  Corrections  are  then  externally  applied.  Experience 
has  shown  that  in  the  dynamic  near-shore  environment,  which  concerns  most 
Corps  projects,  it  is  common  for  soft  sediments  to  be  overlaid  by  more  com¬ 
petent  materials.  Therefore,  it  is  very  important  that  results  be  verified  with 
core  data  and  that  seismic  data  be  collected  concurrently  at  different  frequen¬ 
cies.  Any  future  developments  of  this  method  should  include  dependable 
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algorithms  for  automatically  assessing  the  polarity  of  the  reflection  coefficient 
(see  Chapter  6,  “Recommendations”). 

Assessing  attenuation  across  frequency  band.  Another  approach  to  layer 
identification  involves  detailed  spectrum  analysis  to  detect  possible  frequency 
shifts  in  the  reflection  data  at  each  horizon.  Generally  speaking,  acoustic 
energy  experiences  higher  attenuation  in  more  competent  materials,  such  as 
sands,  than  in  the  softer  silts  and  clays.  The  rate  of  attenuation  in  a  particular 
sediment  unit  is  also  a  function  of  the  frequency,  i.e.,  attenuation  increases 
with  increasing  frequency.  The  frequency  spectrum  from  data  having  good 
S/N  should  indicate  less  energy  at  the  higher  frequencies  for  subbottom  reflec¬ 
tors  below  sand  layers  and  little  to  no  signal  attenuation  across  the  frequency 
spectrum  for  horizons  below  a  soft  sediment  layer.  This  process  has  been  used 
extensively  as  an  aid  in  the  identification  of  sediment  layers.  Figure  39  refers 
to  the  in  situ  condition  of  Figure  38  and  shows  an  example  of  a  clay  pocket 
between  sand  layers.  In  this  case,  the  competent  surface  sands  reflect  most  of 
the  higher  frequency  data.  Little  signal  attenuation  occurs  at  the  bottom  of  the 
clay  layer.  The  center  frequency  of  each  reflection  horizon  is  shown  “shifting” 
to  the  lower  frequencies  below  the  sand  layer,  and  no  comparable  shift  occurs 
below  the  clay  layer. 

To  summarize,  no  one  general  solution  accurately  assesses  all  real-world 
situations.  Without  adequate  core  data,  it  is  important  when  making  predic¬ 
tions  of  sediment  characteristics  that  the  solution  be  derived  from  a  variety  of 
different  methodologies  if  possible.  This  chapter  has  presented  a  robust  suite 
of  tools  that  can  be  successfully  used  to  perform  a  given  characterization 
study.  When  utilized  properly  with  data  of  good  S/N  and  adequate  ground 
truth  information,  the  AI  method  developed  here  can  provide  reliable  infor¬ 
mation  regarding  engineering  properties  of  marine  sediments.  To  illustrate  this 
point.  Chapter  4,  “Case  Studies,”  wiU  present  numerous  examples  of  “ground- 
truthed”  data  from  a  variety  of  field  situations. 
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Input  from  field 
data/tapes 


Figure  2.  Acoustic  impedance  processing  flowchart 
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Figure  3.  Simulated  reflection  energy  envelopes  showing  influence  of  material 
type  on  signal  amplitude 


Figure  4.  Waves  generated  at  an  interface  by  an  incident  P-wave  (after  Officer  (1958)) 
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Atteniiatioa 
l-4b.  Symbols:  drdes-sands 


dau  from  Tables 


imdv^t  sand-^t-dav)*  sand  data  at  500* and  1000  khz  from  Busby  and  Richardson  {mi).  U)w-ircaucncy 
datafline  A-Zemstov,  19^  (land,  sedimentary  strata);  E^d^- 

sand);  C-Bennett,  1967  (sea  floor,  reflection  technique);  Z-a^m  Qii  Vasd ev  Gurevi^  1^2),  E  B 
Epinatyeva  et  al,  and  Berzon  (in  Zemstov,  1969).  Line  labelled/*  indicates  slope  of  any  line  having  a  dependence 
cd  attenuation  on  the  first  power  of  frequency. 


Figure  7.  Attenuation  as  a  function  of  frequency  (Hamilton  1972a) 
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Mean  ^ain  size  mni  versus  Jt  (in  a  =  kp).  Data,  symbols,  and  remarks:  same  as  in  caption  for  figure  3. 
Plot  of  mean  grain  size  in  mm  emphasizes  relationships  in  sands  (see  text).  Regression  equations  for  solid  lines: 
Coarse,  medium,  and  fine  sand,  in  part  (0.6  to  0.167  mm):  i  =0.5374— 0.1 1 13  f.U,).  Fine  sand  (in  part),  and  very 
fine  sand  (0.167  to  0.063  mm) :  i  =  0.S439— 1.9431  (.U,). 


-  .  Forosity  »,  percent,  versus  i  (in  —  Data,  symbols,  and  remarks  same  as  m  caption  for  Figure  3. 

Regression  equations  for  solid  lines:  Coarse,  mediurn,  and  fine  sand  (36  to  46.7  percent):  0.^747+0.00527  (n). 
Very  fine  sand  and  lower-porosity  mixed  sizes  (46.7  to  52  percent):  i  =-0.04903  (n^  1.7688.  Mixed  sizes  (52  to  65 
percent):  i»3J232  -0.0489  (n).  Silt-clays  (65  to  90  percent):  Jk=- 0.7602 -0.0 1487  («)+0.000073  («)*. 


Figure  8.  Attenuation  versus  porosity  and  mean  grain  size  (Hamilton  1972a) 
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ACOUSTIC  IMPEDANCE  ACOUSTIC  IMPEDANCE 


5000 


DELCOAST:  DENSITY  VS  IMPEDANCE 


DELCOAST:  MEAN  GRAIN  SIZE  VS.  IMPEDANCE 


MEAN  GRAIN  SIZE,  PHI 


Figure  10.  Delaware  coast  model:  impedance  versus  density:  impedance  versus  mean  grain 
size 
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Density,  G/Cc 
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10  Channel  Tuned  Circuit  System 


Amplifier  gain  calibration  curve 


TRANSMITTING  RESPONSE 
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Figure  14.  Typical  transducer  specification  data 


Figure  15.  Delaware  coast  survey  geometry  and  computation  of  directivity  index  {Dl)  versus 
water  depth  and  transducer  separation 
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FILE  NO. 
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Figure  16.  Typical  3.5-kHz  caHbration  data 
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Figure  17.  Example  source  level  (SL)  and  bottom  loss  (BL)  calibration 


RECEIVER  SENSITIVITY  CALIBRATION 


Amplifier 


SONAR  EQUATIONS: 

BL  -  SL  -  +  Nhydr  +  ~ 

BL  =  SL  —  Nv/  +  —  SIGq 

COMBINE  AND  SOLVE  FOR  Nhydr 
^  HYOR  “  ^  HYDR  ^  Rr  ~  ^  Rc  ^AR  "T  N  aC 


Figure  18.  Receiving  array  sensitivity  calibration  configuration 
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Calibration  Hydrophone  Signal  FILE:  pci  10043 
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Example  receiving  array  sensitivity  calibration 


Consecutive  Data  FHe  Bl  Versus  Consecutive  Data  Fite 
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Figure  20.  Typical  surface  sediment  characterization 
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Figure  21.  Determination  of  BL  from  multiple  reflections 


Figure  22.  Typical  AC50  output:  amplitude  plot 
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igure  23.  Typical  AC50  output:  impedance  plot 
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Figure  25.  Typical  AC50  output:  alternating  impedance  and  raw  data  plot 


Frequency,  Hz 
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Figure  26.  Typical  spectral  output  (DSA10);  pinger  (3.5-kHz)  reflection  data 
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Figure  27.  Typical  spectral  output  (DSA10):  boomer  reflection  data 


Figure  28.  Typical  spectral  output  (DSA10):  chirped  source  wavelet  data 
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AMPLITUDE 


input  waveform 
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Figure  30.  Correlation  analysis;  boomer  with  noise 
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Figure  31 .  Correlation  analysis:  chirped  signal  with  noise 
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Figure  32.  Correlation  output:  raw  data  versus  correlated  (boomer) 
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DESCRIPTION 


0-29”  Fine  sand,  little  nedlun  sand,  approximately 
15-20%  fine  shell  hash,  increase  in  shell  content 
with  depth 

~  *«<. 

^v.  Wt 

lov*  ct»ii 


.29-32”  Abundant  shell  hash  nixed  with  fine  sand 
32-50”  Mixed  soft  nud,  fine  sand,  some  interbedding 

50-55”  Mixed  fine  sand,  f ine-nedium  ‘shell  hash 

55”  -  bottom  or  core  -  Hard  packed  mud,  dark  grey, 
little  random  mixed  in  very  fine  sand 

Bottom  of  Core 
5'3";  63” 


Table 

Simulated  Site  Core  Data 


Reflection 

Depth 

Density 

Velocity 

Impedance 

Coefficient 

1500 

0.000  • 

>  Water 

3734 

0,409  ■ 

-  Sand 

2198 

-0.258 

'  Soft  mud 

4400 

0,337 

•  Hard  clay 

Sand  over  soft  mud  over  hard  clay 
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I 
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DESCRIPTION 


*  0-25  1/2”  Fine  to  medium  fine  clean  sand,  some 

^.■>4  ’  shell  hash 

to  18”  Band  of  coarse  shell  hash)  « 

4  AA«  -Z'A 

fLV.  *.ll 


’'25  1/2”  -  29  1/2”  Muddy  sand 
■=^~_^'29  1/2"  to  52  1/2"  Very  hard  dewatered  mud 


jr_rz.  Bottom  of  Core 
-  52  1/2" 


Table 

Simulated  Site  Core  Data 


(Reflection 
Coefficient 


1.00 

1500 

2.03 

1836 

2.20 

2000 

Sand  over  hard  clay 


File  Reflno3.asc  Reflection  Amplitude  vs.  Depth  -  ChpSmd 
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Figure  34.  Synthetic  seismograms:  sand/soft  mud/hard  clay 


Synthetic  seismograms:  sand/hard  clay 


Matched  Filter  Correlation  -  DPC10  Matched  Filter  Correlation  -  DPC10 


Figure  36.  Synthetic  seismogram  matched  filter  correlation  -  boomer  and  chirper 
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Reflection  Coefficient 


GEOACOUSTIC  MODEL  -  LINE  7 
ACOUSTIC  PROPERTIES 
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Note:  Seismogram  in  depth  (ft),  boomer  data  in  travel  time. 


Figure  38.  Forward  model  comparison  with  actual  data 
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Communication 


Technology  transfer  as  illustrated  by  Ballard  and  Grau  (1992)  and  Sjostrom, 
Ballard,  and  McGee  (1991)  in  brochures  and  videos  shows  how  AI  surveys 
can  help  define  bottom  and  subbottom  conditions  prior  to  dredging.  Once  a 
need  is  established,  planning  unfolds.  The  necessity  for  two-way  communi¬ 
cation  between  tlie  project  manager  and  the  scientist/engineer  conducting  an 
acoustic  impedance  survey  cannot  be  overemphasized.  For  maximum  benefits 
to  be  derived,  this  should  take  place  at  a  point  early  in  the  planning  stages  of 
a  project,  rather  than  in  the  latter  stages  after  a  boring  exercise  is  already 
underway.  A  face-to-face  meeting  between  all  involved  parties  at  the  project 
office  should  take  place,  as  opposed  to  solely  written  or  telephone  communica¬ 
tion.  A  site  visit,  if  practical,  is  desirable.  In  the  course  of  the  onsite  meeting, 
the  project’s  background  is  presented  by  the  project  manager  along  with  both 
primary  and  secondary  objectives.  Potential  problem  areas  must  be  accounted 
for  and  limitations  of  survey  results  addressed.  After  objectives  are  agreed 
upon,  site  layout  and  optimum  survey  grids  can  (and  must)  be  established  prior 
to  making  estimates  of  time  and  cost.  The  time  required  to  perfomi  the  survey 
is  directly  dependent  on  the  number  of  miles  to  be  covered  and  operational 
factors  such  as  weather,  shipping  traffic,  and  other  site-dependent  obstacles, 
and  the  survey  cost  is  dependent  on  time  and  survey  vessel  support. 

The  single  most  elusive  factor  influencing  time  and  cost  is  weather.  Oper¬ 
ating  conditions  must  be  safe  and  conducive  to  gathering  high  signal-to-noise 
ratio  data.  In  these  terms,  “sea  state”  is  governed  by  such  factors  as  where 
operations  are  conducted— in  open  water  as  opposed  to  a  sheltered  harbor,  in 
deep  (>100  ft)  as  opposed  to  shallow  (<30  ft)  water.  In  general,  wave  height 
should  not  exceed  about  3  ft  for  acceptable  operations  because  pitch  and  roll 
of  the  transducers  and  survey  platfonn  can  produce  unknown  signal  direction¬ 
ality.  Whitecapping,  which  is  a  function  of  wind,  water  depth,  and  shelter 
conditions,  also  produces  excessive  seismic  noise.  For  these  reasons,  opera¬ 
tional  decisions  must  be  made  (from  Uie  safety  standpoint)  on  a  daily  basis  by 
the  survey  vessel  pilot  in  conjunction  with  the  data  acquisition  expert. 

Weather  allowances  are  primarily  based  upon  seasonal  history  but,  in  the  view 
of  both  parties,  should  be  subject  to  change.  In  the  early  stages  ol 
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communication  and  negotiation,  agreement  should  be  reached  concerning  the 
handling  of  standby  or  “weather”  days. 

A  comprehensive  review  of  background  materials  including  existing  boring 
logs  and  site  geology  must  also  be  performed.  Support  requirements  of  the 
survey  team  must  be  addressed  along  with  the  time  frame  for  survey  perfor¬ 
mance,  interim  and  final  reporting,  and  costs. 

All  of  the  above  factors  are  critical  to  optimum  job  performance  and  bene¬ 
fits.  To  illustrate  the  manner  in  which  good  survey  planning  should  occur,  an 
actual  example  of  dialogue  exchange  between  the  New  England  Division 
(NED)  and  WES  is  presented  in  Appendix  B.  This  project,  fully  documented 
by  Sjostrom  and  Leist  (1994),  was  chosen  not  only  tecause  of  its  overall  suc¬ 
cess,  but  as  an  example  of  support  provided  by  the  Dredging  Operations  Tech¬ 
nical  Support  (DOTS)  Program.  Documents  used  for  illustration  include  an 
information  summary  sheet  (Figure  Bl)  on  the  Boston  Harbor,  MA,  Navigation 
Improvement  Project  which  provided  the  necessary  background  to  WES  per¬ 
sonnel  for  preparation  of  a  face-to-face  plaiming  meeting  and  seminar  held  at 
NED.  Also  included  is  a  Memorandum  for  Record  (Figure  B2)  prepared  by 
WES  scientists  documenting  that  meeting.  The  final  proposal  submitted  by 
WES  to  NED  is  also  included  (Figure  B3)  because  of  its  detailed  treatment  of 
such  items  as  the  following: 

a.  A  clear  statement  of  the  objective. 

b.  The  scope  of  work  to  be  performed  (including  instrumentation). 

c.  The  proposed  survey  grid  to  be  used  at  four  sites  of  interest. 

d.  Support  requirements  to  be  provided  by  NED. 

e.  Scheduling  and  reporting  (both  interim  and  final). 

These  serve  to  further  illustrate  the  point  that  communication  is  critical 
throughout  the  entire  planning  stage. 


Surveying  and  Reporting 

During  field  experience  with  conduct  of  AI  surveys  some  pitfalls  became 
readily  apparent.  Foremost  were  navigational  problems.  Obviously,  data  are 
worthless  if  locations  are  not  known;  consequently,  great  care  must  be  taken  to 
ensure  that  accurate  grid  coordinates  are  finnly  tied  to  the  seismic  data  acquisi¬ 
tion  system,  i.e.,  the  systems  must  be  able  to  communicate.  Early  problems 
arose  as  the  result  of  interface  compatibility  between  navigation  equipment 
provided  by  the  support  element  and  the  seismic  data  acquisition  system. 

Once  recognized,  these  problems  were  quickly  overcome  but  required  a 
thorough  understanding  of  the  system  to  be  used-its  setup  time  requirements 
(range-range  land  stations,  manpower,  etc.)  and  limitations  (in  some  instances. 
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data  were  not  stored  by  survey  instrumentation,  and  if  problems  arose,  an 
entire  seismic  survey  line  would  have  to  be  repeated).  A  second  pitfall  was 
establishment  of  the  water  depth  datum  plane.  Many  Corps  support  elements 
use  200-kHz  echosounders  exclusively  to  define  “bottom.”  In  actuality  this 
can  be  the  top  of  a  fluid  mud  layer  rather  than  true  hard  bottom  as  defined  by 
the  AI  seismic  source.  A  solution  was  reached  when  echo  sounder  data  were 
compared  to  pinger  data  and  the  fluid  mud  thickness  accounted  for.  In  final 
analysis,  the  actual  datum  plane  or  bottom  elevation  is  established  with  refer¬ 
ence  to  tide  gauges  in  accordance  with  local  operational  procedures.  No 
longer  does  either  of  these  pitfalls  involving  X,Y,Z  coordinates  present  a  prob¬ 
lem.  The  WES  AI  integrated  data  management  and  navigation  system  coor¬ 
dinates  satellite  (GPS)  navigation  with  an  Innerspace  445  echo  sounder  and 
seismic  source  data.  All  combined  data  are  now  stored  on  hard  disk. 

Visualization  of  critical  data  will  be  an  important  factor  for  the  user,  whose 
preferences  should  be  made  explicit  as  exemplified  by  NED.  In  that  specific 
case,  the  grid  shape  would  not  have  been  conducive  to  three-dimensional  (3-D) 
visualization,  and  users  actually  preferred  cross-sectional  representations  of 
channel  conditions. 

Actual  conduct  of  the  survey  must  be  planned  well  in  advance  so  that  data 
can  be  acquired  in  a  manner  consistent  with  the  desired  type  of  presentation. 
Data  presentations  most  often  consist  of  two-dimensional  (2-D)  cross  sections 
color-coded  in  terms  of  material  densities  or  type  of  material.  (In  earlier  pre¬ 
sentations,  2-D  plots  in  shades  of  gray  were  the  only  choice.)  However,  the 
user  may  wish  to  color-code  a  site  in  three  dimensions  so  that  “pockets”  of 
materials  and  volumetric  calculations  take  on  additional  meaning.  It  should  be 
noted  that  3-D  presentations  are  more  time-consuming  and  consequently  more 
costly  to  produce,  but  in  the  words  of  many  users  “well  worth  the  price.” 
Examples  of  various  types  of  visualization  will  be  discussed  and  illustrated  in 
Chapter  4. 

If  at  all  possible,  the  project  manager  or  the  person(s)  using  the  survey  data 
should  physically  participate  in  the  actual  field  exercise  and  data  processing/ 
interpretation  phases  of  the  survey,  thereby  becoming  an  integral  part  of  the 
“team”  and  developing  a  keen  insight  into  the  quality  of  data  and  its  limita¬ 
tions.  Further,  many  decisions  can  and  should  be  made  onboard  the  survey 
vessel  as  data  are  acquired. 

Upon  completion  of  the  final  report,  a  comprehensive  discussion  of  findings 
should  be  presented  in  person  by  the  surveyors  to  the  project  manager  and 
other  end-product  users  in  their  office.  This  discussion  should  include  all 
aspects  of  the  survey  from  planning  through  conduct,  data  processing,  inter¬ 
pretation,  computer  biases,  and  limitations.  Recommendations  should  be  made 
regarding  placement  of  verification  borings  and  the  degree  to  which  objectives 
have  been  met. 
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4  Case  Studies 


Overview 


Presented  in  this  section  are  site  calibration  and  verification  results  (ground 
truthing)  conducted  at  a  number  of  test  sites  around  the  country.  Although 
surveys  have  been  conducted  at  some  16  sites,  ground  truth  information  is  not 
available  for  all  sites.  In  some  instances,  sampling  is  still  underway.  Hence, 
only  those  data  obtained  from  projects  whose  objectives  fully  exploited  the  AI 
procedures  wiH  be  presented.  Only  two  examples,  the  Gulfjxrrt  fluid  mud 
study  and  one  area  of  Oakland  Harbor,  were  actually  funded  by  DRP.  Other 
surveys,  having  a  variety  of  objectives,  were  performed  on  a  fully  reimbursable 
basis,  sponsored  by  various  Corps  District  and  Division  offices  or  other 
agencies.  Geoacoustic  data  are  presented  for  direct  comparison  to  ground  truth 
information  derived  from  conventionally  accepted  Corps  practices.  These 
comparisons  consist  of  acoustically  derived  sediment  properties,  and  measured 
sediment  properties  and/or  lithology.  Since  it  was  the  ultimate  goal  of  this 
work  to  produce  a  validated  procedure  sufficiently  accurate  for  characterizing 
bottom  and  subbottom  materials,  a  concerted  effort  was  made  to  quantify  the 
degree  of  accuracy.  It  was  recognized  that  the  procedures  should  provide 
unrestricted  use  by  Corps  personnel  for  planning  purposes  or  even  as  evidence 
for  use  in  court  dredging  claims.  In  virtually  all  instances,  geoacoustic  predic¬ 
tions  were  accurate  to  within  ±10  percent  while  in  many  cases  accuracy  was 
within  ±  5  percent. 

In  addition  to  discussions  on  ground  truth,  examples  of  2-  and  3-D  visual¬ 
ization  and  volumetric  calculations  are  presented.  These  serve  to  provide 
potential  users  with  the  versatility  of  this  technique. 

Following  is  a  complete  listing  of  sites  surveyed  to  date  using  the  acoustic 
impedance  approach: 
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Project 

Sponsor 

Oakland  Harbor 

San  Francisco  District 

Mobile  Ship  Channel 

Mobile  District 

Gulljxrrt  Ship  Charmel 

Mobile  District 

Panama  City  Coastline 

Mobile  District 

Pascagoula  Ship  Charmel 

Mobile  District 

Savannah  Harbor 

Savannah  District 

Galveston  Ship  Qiannel 

Galveston  District 

Houston  Ship  Charmel 

Galveston  District 

Delaware  River  Charmel 

Philadelphia  District 

Delaware  Bay  Coastlines 

Philadelphia  District 

New  Jersey  Shoreline 

Philadelphia  District 

Delaware  River  Ship 

Charmel/  Harbor 

Philadelphia  District 

Boston  Harbor 

New  England  Division 

Baltimore  Harbor 

Baltimore  District 

Picatiimy  Lake 

Army  Environmental  Center 

Clinch  River/Poplar  Creek 

Department  of  Energy 

Charleston  Harbor/Ship  Charmel 

Charleston  District 

For  reasons  previously  explained,  it  is  beyond  the  scope  of  this  report  to 
present  all  data;  however,  substantive  data  from  a  wide  range  of  geologically 
different  sites  will  be  discussed  and  statistically  reported. 


Gulfport  and  Mobile  Ship  Channels 

The  Gulfix)rt  field  exercise  was  conducted  in  the  summer  of  1990  as  the 
second  evaluation  of  the  acoustic  impedance  method  developed  under  DRP 
sponsorship.  During  this  evaluation  phase,  the  target  application  was  to 
explore  the  effectiveness  of  the  technique  in  density  determinations  of  fluid 
mud.  Figure  40  shows  amplitude  plots  of  20-kHz  and  45-kHz  surveys  con¬ 
ducted  in  Gulfport  Ship  Charmel,  MS.  Figure  41  shows  acoustic  impedance 
density  predictions  versus  nuclear  densitometer  data  for  fluff/fluid  mud  type 
materials  in  the  Gulfport  Ship  Channel  as  shown  in  Ballard,  McGee,  and 
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Whalin  (1992).  Correlation  proved  to  be  excellent,  within  an  average  percent¬ 
age  difference  of  -1.60. 

Data  acquisition.  Under  a  reimbursable  agreement  with  Mobile  District, 
Gulfix)rt  Ship  Channel  was  revisited  to  survey  the  route  of  a  proposed  new 
channel.  The  primary  objective  of  this  investigation  was  to  determine  densities 
of  subbottom  materials,  isolate  sands,  and  compute  an  estimate  of  their  vol¬ 
ume.  Mobile  Ship  Channel  was  investigated  separate  from  Gulfport,  but  since 
the  projects  are  geographically  similar  and  were  both  sponsored  by  the  Mobile 
District,  it  is  fitting  to  describe  their  results  collectively.  However,  prior  to  a 
discussion  of  the  ground  truthing  of  these  data,  a  comparison  with  Hamilton’s 
worldwide  database  will  prove  interesting.  A  plot  of  the  impedance  function 
versus  laboratory  measurements  of  density  from  core  samples  taken  in  the 
Mississippi  Sound  and  Mobile  Bay  is  presented  in  Figure  42.  A  slight 
modification  by  Caulfield  (1992a)  to  Hamilton’s  data  is  represented  by  the 
solid  line  in  Figure  42.  It  shows  remaiicable  agreement  with  the  data  at  both 
sites.  Figure  43  shows  the  acoustic  predictions  versus  core  data  for  consoli¬ 
dated  materials  in  Mobile  Ship  Channel.  In  this  case,  results  are  within  1  per¬ 
cent.  Figure  44  shows  the  same  data  for  materials  in  Gulfport  Ship  Channel. 
Note  that  this  comparison  format  (Figure  44)  is  currently  preferred  by  users. 
This  figure  displays  a  segment  of  the  actual  acoustic  reflection  record 
color-coded  to  represent  changes  in  signal  amplitude  followed  by  the  acoustic 
impedance  processing  of  that  record.  Next,  the  computed  density  plotted  as  a 
function  of  depth  is  shown  along  with  laboratory-derived  densities  from  core 
samples  in  the  same  vicinity.  Finally,  the  predicted  lithology  computed  from 
acoustic  data  is  shown  adjacent  to  lithology  obtained  from  core  analysis.  (If 
no  core  samples  are  available,  the  core  section  is  omitted  from  the  display.) 

Visualization.  Three  2-D  profile  views  representative  of  a  segment  of  the 
Gulfport  Ship  Channel  are  illustrated  in  Figure  45  with  survey  line  9  posi¬ 
tioned  along  the  center  line  of  the  proposed  channel.  Lines  1  and  18  are  paral¬ 
lel  to  each  edge  of  the  proposed  channel.  This  type  of  display  delineates  the 
extent  of  pertinent  density  zones  of  interest  to  the  engineer,  and  the  virtually 
continuous  data  coverage  greatly  decreases  the  possibility  of  significant  mate¬ 
rial  changes  going  undetected. 

Density  predictions.  By  incorporating  the  virtually  continuous  coverage  of 
subbottom  materials  with  digital  terrain  modeling  techniques,  rapid  and  accu¬ 
rate  computations  can  be  made  of  volume  and  material  type  to  be  removed  by 
dredging.  Furthermore,  a  detailed  database  has  now  been  established  for  proj¬ 
ect  monitoring  and  long-term  planning.  Computed  sediment  densities  within 
the  project  area  can  be  displayed  in  a  color-coded,  3-D  view  as  shown  in  Fig¬ 
ure  46  if  desired  by  the  user.  In  this  example,  red,  brown,  and  yellow  shad¬ 
ings  are  indicative  of  less-dense  material.  Blues  represent  sands. 

Volumetric  calculations.  At  the  project  planner’s  discretion,  an  area  of 
interest  may  be  viewed  from  other  angles  or  other  displays  created  by  stripping 
or  slicing  at  any  desired  coordinate.  Volume  of  any  material  to  be  removed 
can  be  easily  calculated  and  displayed  as  shown  in  Figure  47.  This  example 
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predicts  the  configuration  of  sands  underlying  the  segment  of  Gulfport  Ship 
Channel  shown  in  Figure  46.  Ground  truth  information  concerning  the  Gulf¬ 
port  Ship  Channel,  as  related  to  predictions  of  material  type  and  location, 
could  only  be  obtained  qualitatively  from  the  Mobile  District  project  manager 
as  dredging  progressed.  In  his  words,  the  materials  and  quantities  predicted 
closely  matched  those  dredged. 

It  should  also  be  noted  that,  possibly  for  the  first  time  in  USACE  history. 
Mobile  District  advertised  site  conditions  for  dredging  at  Gulfport  Ship  Chan¬ 
nel  using  density  predictions  determined  by  geophysical  acoustic  impedance 
surveying  in  addition  to  information  normally  supplied.  Infonnation  derived 
from  continuous  site  characterization  rather  than  borings  alone  is  considered 
(by  Mobile  District  personnel)  to  be  beneficial  to  the  potential  dredging  con¬ 
tractor  and  the  Corps  for  planning  purposes. 


Galveston  Ship  Channel 

During  the  summer  of  1992,  an  acoustic  impedance  survey  of  Galveston 
Ship  Channel  was  conducted  (synopsized  in  Ballard  et  al.  (1992)).  Although 
several  objectives  were  accomplished,  the  phase  concerning  volumetric  deter¬ 
minations  of  “dirty”  sands  (density  of  1.7-1.95  g/cc)  underlying  the  proposed 
channel  was  of  specific  interest  because  of  its  potential  use  for  beaches.  The 
following  discussion  presents  the  methodology  used  to  locate,  identify,  and 
quantify  that  material  of  interest  to  the  U.S.  Amiy  Engineer  District, 

Galveston. 

Data  acquisition.  Figure  48  is  a  typical  3.5-kHz  pinger  amplitude  cross 
section  reproduced  from  a  color  field  record  obtained  in  Galveston  Ship  Chan¬ 
nel.  Note  that  the  top  of  the  graph  is  not  the  water  surface  but  an  assigned 
water  column  delay.  This  offset  allows  full  vertical  expansion  of  the  sub¬ 
bottom  display,  which  in  this  case  extends  into  the  subbottom  more  than  40  ft. 
Changes  in  stratigraphy  are  readily  apparent.  Records  of  this  type  are  used  as 
quick-look  guidance  in  boring  placement. 

Density  predictions.  Figures  49-53  display  acoustic  predictions  versus 
core  data  for  consolidated  materials  in  Galveston  Ship  Channel.  Analysis 
shows  acoustically  derived  density  values  vary  +0.2  g/cm^  from  in  situ  values. 
As  previously  described,  computed  sediment  densities  within  the  project  area 
can  be  displayed  in  a  color-coded,  3-D  view  (Figure  54).  In  this  example, 
lighter  shadings  are  indicative  of  less-dense  material-the  darkest  analogous  to 
densities  >1.7  g/cc. 

Volumetric  calculations.  Using  Figure  54  as  an  example,  the  volume  of 
any  material  to  be  removed  can  be  easily  calculated.  Tlris  example  wiU  predict 
the  configuration  of  silty  (“dirty”)  sands  (>1.7  g/cc  density;  2864-3052  x 
10^  g/cm^  sec  acoustic  impedance)  underlying  a  selected  segment  of  the 
project  study  area.  The  volume  of  material  present  within  a  selected  area  of 
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the  perspective  model  is  arrived  at  by  calculating  the  volume  of  material  pre¬ 
sent  within  the  corresponding  area  of  each  profile  line  model. 

In  the  example,  the  profile  line  model  shown  in  Figure  55  corresponds  to 
the  section  of  the  3-D  perspective  model  (see  Figure  54)  nearest  the  viewer. 
Note  that  the  axis  labels  in  Figure  54  show  northing  and  easting  coordinates, 
while  the  axis  label  in  Rgure  55  reflects  relative  feet  from  start-of-line.  This 
conversion  is  necessary  in  order  to  correlate  the  profile  line  models  to  ampli¬ 
tude  records  which  are  plotted  with  locations  relative  to  start-of-line.  Before 
volumes  are  calculated,  the  area  of  interest  must  be  sliced  out  of  the  computer 
model  and  the  material  density  range  to  be  displayed  must  be  selected.  In  this 
example,  the  volume  of  material  present  in  the  area  between  the  bottom  to 
-55.0-fl  depth  and  between  3,000  to  7,000  ft  from  start-of-line,  having  a  den¬ 
sity  of  >1.7  g/cc,  will  be  selected.  The  section  of  the  model  below  -55.0  ft  is 
stripped  away,  and  the  sections  from  0  to  3,000  ft  and  7,000  to  19,000  ft  are 
sliced  from  the  display  as  shown  in  Figure  56.  Finally,  the  density  of  material 
to  be  displayed  is  set  to  >1.7  g/cc.  Calculations  based  on  these  parameters  are 
then  computer-generated  and  displayed  as  shown  in  Figure  57,  yielding  an  esti¬ 
mated  volume  of  167,366  cu  yd  of  “dirty”  sand  (materials  >1.7  g/cc)  within 
the  specified  boundaries.  This  step  is  repeated  for  each  profile  line  model 
within  the  area  of  interest,  summing  the  estimated  volumes  for  an  overall 
project  volume  total. 


Boston  Harbor 

At  the  request  of  NED,  WES  perfonned  a  reimbursable  waterborne  seismic 
reflection  and  side-scan  sonar  survey  of  three  tributaries  to  Boston  Harbor, 
Massachusetts,  during  November  1992.  The  survey,  reported  in  its  entirety  by 
Sjostrom  and  Leist  (1994)  and  further  discussed  in  Appendix  B,  included 
Mystic  and  Chelsea  Rivers,  Reserved  Channel,  and  the  Inner  Confluence  Area. 
The  objective  of  the  survey  was  to  quantify,  as  a  function  of  depth,  bottom  and 
subbottom  sediments  in  terms  of  density  and  soil  type  to  elevation  -42  ft  mean 
low  water  (mlw)  in  all  areas  except  the  Chelsea  River  where  the  desired  datum 
was  -40  ft  mlw.  Information  was  to  be  used  to  facilitate  positioning  of  soil 
borings,  particularly  in  areas  of  suspected  rock  outcroppings.  During  earlier 
phases  of  maintenance  dredging  woric,  NED  had  obtained  some  boring  infor¬ 
mation,  but  little  was  documented  regarding  laboratory  density  analysis.  An 
additional  six  shallow  penetration  cores  were  collected  following  the  geophysi¬ 
cal  survey  to  more  accurately  calibrate  the  geoacoustic  model  and  develop  the 
necessary  acoustic  parameters  to  derive  estimates  of  bottom  and  subbottom 
material  density.  Laboratory  densities  were  obtained  from  these  cores.  Length 
of  the  core  samples  ranged  from  1.6  to  2.0  ft.  Consequently,  direct  com¬ 
parisons  with  acoustic  data  were  restricted  to  the  upper  2.0  ft.  Fifteen  addi¬ 
tional  cores  (containing  no  laboratory  density  infonnation)  ranging  between  6 
and  17  ft  were  completed  in  June  1993  and  were  also  used  for  subbottom 
material  comparison  purposes. 
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Laboratory-derived  density  values  ranged  from  1.38  to  1.84  g/cm^.  Fig¬ 
ure  58  is  a  typical  calibration  comparison  display  showing  both  acoustic  and 
shallow  laboratory-derived  density  values.  This  figure  was  selected  for  presen¬ 
tation  because  it  also  contained  subbottom  lithographic  information  described 
in  NED  boring  FD-93-K.  Note  the  agreement  between  lithology  and  acoustic 
densities,  i.e.,  high  densities  imply  GPs  and  SPs.  Acoustically  derived  density 
estimates  using  the  determined  calibration  parameters  showed  excellent  correla¬ 
tion  with  one  possible  exception  which  showed  the  acoustic  data  to  be  approxi¬ 
mately  10  percent  lower  than  the  laboratory.  In  the  other  cases  a  match  of 
about  5  percent  was  observed.  It  was  interesting  to  note  that,  in  all  cases, 
laboratory  densities  in  the  subbottom  upper  2.0  ft  were  higher  than  those 
obtained  acoustically  in  the  Boston  area. 


Delaware  Atlantic  Coast 


An  acoustic  impedance  survey  sponsored  by  the  Philadelphia  District  was 
conducted  along  the  Delaware  Atlantic  coast  in  the  fall  of  1992  as  part  of  a 
feasibility  study  for  shore  protection  solutions.  This  survey  was  selected  for 
discussion  because  its  objective  was  quite  different  from  that  of  typical  harbor 
dredging  projects.  The  objective  was  to  quantify,  versus  depth,  bottom  and 
subbottom  sediments  in  terms  of  in  situ  density,  mean  grain  size,  and  soil  type 
to  a  depth  of  about  20  ft  below  the  ocean  bottom,  where  possible,  for  a 
3-mile-wide  area  between  Cape  Henlopen  and  Fenwick  Island  along  the  Dela¬ 
ware  coast.  The  results  were  to  provide  initial  estimates  of  the  sediment  char¬ 
acteristics  for  the  purpose  of  defining  the  limits  of  available  granular  materials 
which  could  be  used  for  a  variety  of  purposes. 

Geoacoustic  calibration.  The  geoacoustic  calibration  was  accomplished 
using  sediment  grain  size  versus  acoustic  impedance.  Density  values  were 
computed  based  on  empirical  relationships  for  density  and  mean  grain  size, 
such  as  presented  in  Figure  11.  Figure  59  presents  a  summary  of  the  geoacou¬ 
stic  calibration  of  3,500-Hz  pinger  data  with  core  information  collected  for  the 
Delaware  coast  survey.  Each  figure  includes  a  representative  section  of  the 
reflection  data  in  the  vicinity  of  the  core;  plots  of  the  acoustically  derived 
estimates  of  impedance,  density,  and  mean  grain  size  versus  depth;  measured 
mean  grain  size  versus  depth;  and  the  core  lithology.  Note  that  the  format  of 
these  plots  differs  from  the  color  versions  presented  previously,  illustrating  the 
fact  that  visual  data  presentation  must  be  tailored  to  the  needs  of  the  user.  In 
this  case,  shades  of  gray  suffice  to  convey  the  message  intended.  Density 
predictions,  based  upon  Hamilton’s  extensive  database,  can  only  be  estimated 
to  be  within  ±10  percent  of  in  situ  since  no  density  measurements  were  con¬ 
ducted  on  the  core  samples.  According  to  an  assessment  by  Sjostrom  and 
Leist  (1994),  had  laboratory  density  analysis  been  available  for  calibration 
purposes,  these  results  could  likely  have  been  improved  to  possibly  5  to  7  per¬ 
cent.  Estimates  of  mean  grain  size,  especially  for  the  surface  sediments,  are 
within  about  ±  0.54),  based  upon  actual  laboratory  analysis.  This  is  quite 
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encouraging  since  acoustic  response  characteristics  of  sediments  as  related  to 
mean  grain  size  tend  to  exhibit  more  scatter  in  the  data. 

Example  sediment  profile.  After  calibration  and  verification,  the  acoustic 
reflection  data  are  processed  to  characterize  the  sediment  properties  throughout 
the  survey  area.  An  example  presentation  of  the  distributions  of  computed 
sediment  densities  and  mean  grain  sizes  within  the  project  area  is  given  in  Fig¬ 
ure  60  as  a  2-D  profile  illustrating  the  primary  bottom  and  subbottom  inter¬ 
faces  and  differing  zones  of  sediment  material.  The  profile  illustrates  the 
depth  to  a  particular  interface  in  feet  mean  lower  low  water  (mllw),  represen¬ 
tative  sediment  properties,  and  corresponding  location  along  the  survey  line. 

An  area  of  primarily  granular  sediments  is  identified  near  the  southern  end  of 
the  survey  line.  Two  relic-incised  channels  filled  with  lagoonal  muds  and  silts, 
and  reworked  clayey  sands,  are  detected.  The  labelled  black  dots  at  the  top  of 
each  profile  denote  the  survey  track  line  and  direction.  Each  dot  also  repre¬ 
sente  the  beginning  of  every  seismic  data  file  recorded  in  order  to  give  an 
indication  of  the  data  coverage  along  each  line  and  assist  in  correlating  the  raw 
data  and  interpreted  results.  The  associated  label  represents  the  data  file  num¬ 
ber  and  correlates  with  the  data  file  number  on  the  color-coded  subbottom 
reflection  records.  Lines  of  latitude  are  displayed  on  each  profile.  The  sedi¬ 
ment  profiles  have  been  completely  adjusted  for  horizontal  position  (effects  of 
boat  speed)  and  survey  heading.  All  profiles  are  presented  heading  in  a  south¬ 
erly  direction,  allowing  consistency  in  the  data  interpretation.  Actual  boat 
heading  is  in  the  direction  of  increasing  data  file  number?  on  the  profiles. 
Finally,  all  cores  collected  along  the  survey  line  are  identified  and  core  logs 
plotted  directly  on  the  profiles. 


Oakland  Harbor 


An  abbreviated  AI  survey  of  Oakland  Harbor  was  conducted  in  conjunction 
with  a  conventional  subbottom  survey  which  was  ongoing  as  a  reimbursable 
project  by  WES  for  the  San  Francisco  District.  The  objective  of  the  survey 
was  to  differentiate  “young”  bay  muds  from  more  competent  sand  sediments 
for  expeditious  environmental  testing.  Circumstances  related  to  the  ongoing 
survey  provided  a  “target  of  opportunity”  to  simultaneously  collect  data  in  the 
AI  format.  Limited  DRP  funds  were  allocated  for  this  study  to  evaluate  the  AI 
concept  in  a  highly  complex  geologic  environment.  Figures  61-64  show  a 
comparison  of  the  geoacoustic  data  and  boring  lithology.  Of  particular  interest 
on  most  of  the  density  plots  (which  continue  to  evolve  in  format  according  to 
sponsor  request)  is  a  plot  labeled  “sign  acoustic  core.”  These  data  are  output 
from  a  newly  developed  software  program  designed  to  automatically  determine 
the  polarity  of  the  reflection  coefficient.  It  is  included  for  comparison  with  the 
“standard”  method  of  acoustic  core  determinations.  This  new  method  for 
assessing  correct  polarity  appears  to  work  reasonably  well,  but  comparisons 
will  continue  to  be  monitored  until  sufficient  detailed  analytical  verification 
has  been  accomplished.  Figure  65  presents  a  typical  sediment  profile  from  the 
Oakland  Harbor  Iimer  Tangent.  The  sediment  profile  developed  using  the  AI 
method  displays  a  layered  sediment  environment  consisting  of  muds,  clays,  and 
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silts  overlaying  the  more  competent  sands.  Figure  66  is  a  color-coded  (by 
material  type)  3-D  visualization  of  the  surveyed  area. 


Statistical  Assessment 

Figure  67  shows  the  correlation  of  laboratory-measured  densities  from  core 
samples  with  acoustically  derived  estimates  of  density  at  those  locations.  It 
must  be  emphasized  that  most  of  the  samples  are  disturbed  samples  collected 
using  vibracore  methods  which  can  result  in  possible  errors  with  the  laboratory 
measurements.  Even  so,  the  95 -percent  confidence  interval  (i.e.,  95  percent  of 
the  data  fall  within  this  boundary)  is  ±  0.19  g/cm^.  The  standard  error  esti¬ 
mate  of  acoustic  density  predictions  as  compared  to  the  “perfect  agreement”  is 
±  0.10  g/cm^  .  In  other  words,  data  verify  the  fact  that  density  can  be  acousti¬ 
cally  derived  to  within  ±  10  percent  of  in  situ  conditions. 
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Figure  40.  Amplitude  color  plot,  20-  and  45-kHz  survey  profiles 


Density  Estimation  Comparison 
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Figure  41 .  Acoustic  impedance  density  prediction  versus  nuclear  densitometer 
data  for  fluid  mud 
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Figure  42.  Hamilton’s  impedance  function  (modified  by  Caulfield  (1992a))  with 
Gulfport  and  Mobile  data  superimposed 
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Density  Estimations  -  Mobile  Bay 
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Figure  43.  Acoustic  predictions  versus  core  data  -  Mobile 
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Rgure  44.  Acoustic  predictions  versus  core  data  -  Gulfport  (Sheet  1  of  5) 
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Figure  44.  (Sheet  4  of  5) 
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Figure  45.  Density  cross  sections  in  Gulfport  Ship  Channel 
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Figure  46.  3-D  view  of  a  section  of  Gulfport  Ship  Channel 


Figure  47.  Volumetric  calculations  for  sand  in  a  section  of  Gulfport  Ship  Channel 
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LOCATION:  Core  CERC-4,  Anchorage  Basin,  Galveston,  TX 
POSITION:  WES  Survey  Line  AP08,  approx,  Station  17+000 
EASTING:  3359335 
NORTHING:  576465 
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Figure  49.  Bottom  surface  calibration  at  WES  Core  CERC-4,  Anchorage  Basin 
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LOCATION:  Cores  CERC-3  <Sc  3ST-87,  Anchorage  Basin,  Galveston,  TX 
POSITION:  WES  Survey  AP04,  approx.  Station  16+500 
EASTING:  3359637 
NORTHING:  575490 
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Figure  51.  Subbottom  sediment  calibration  at  Cores  CERC-3  and  3ST-87,  Anchorage  Basin 
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Figure  52.  Check  calibration  at  Core  CERC-9,  Inner  Bar  Channel 
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Depth 


Figure  56.  2-D  profile  showing  area  of  interest  within  depth  and  location  boundaries 


Figure  57.  2-D  profile  model  showing  only  material  >1.7  g/cm^  density  within  area  of  interest 
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Figure  58.  Typical  calibration  comparison  display  of  acoustic  data  versus  laboratory  density 
and  lithographies  (Boston  Harbor) 
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GEDACDUSTIC  CALIBRATION 
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Figure  59.  (Sheet  3  of  5) 
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Figure  59.  (Sheet  5  of  5) 


Figure  60.  2-D  profile  of  sediment  material  (Delaware  coast) 
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Figure  61 .  Oakland  Harbor  geoacoustic  calibration  (CORE  CH-3) 
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Figure  62.  Oakland  Harbor  geoacoustic  calibration  (CORE  IC-18) 
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Rgure  66.  3-D  visualization  of  a  portion  of  Oakland  Harbor 
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Figure  67.  Statistical  analysis  of  measured  density  versus  acoustically  derived  density 
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5  Boundary  Conditions  and 
Limitations 


The  method  described  herein  represents  an  engineering  solution  to  the  prob¬ 
lem  of  remotely  assessing  the  physical  characteristics  of  marine  sediments. 

The  system  is  not  a  “black  box”  device  capable  of  assessing  every  conceivable 
geoacoustic  situation  occurring  in  the  real  world;  therefore,  it  is  important  that 
boundary  conditions  and  limitations  of  the  technique  be  emphasized. 


Boundary  Conditions 

A  number  of  important  assumptions  and  limiting  conditions  have  been 
required  to  expedite  a  solution.  These  boundary  conditions  and  assumptions 
employed  in  the  Acoustic  Impedance  method,  noted  throughout  this  report,  are 
summarized  as  follows; 

a.  Sound  wavefront  propagated  as  a  plane  wave  acting  at  normal  inci¬ 
dence  to  a  horizontally  layered  system.  Virtually  all  commercially 
available  imderwater  profiling  systems  produce  spherical  wave  fronts 
with  either  highly  directional  or  omnidirectional  beam  patterns.  Certain 
beam  pattern  and  bottom  type  combinations  can  significantly  influence 
the  quality  of  a  reflected  signal.  For  Coips  projects  such  as  ship  chan¬ 
nels,  disposal  mounds,  or  high-energy  sediment  zones,  the  bottom 
topography  may  not  be  horizontal.  A  rapidly  changing  bottom  topog¬ 
raphy  will  alter  the  integrity  of  returned  echoes  by  either  redirecting  the 
echo  away  from  the  receiver  or  focusing  too  much  energy  towards  the 
receiver.  Also,  side  echoes  reflected  off  vertical  objects  or  barriers  can 
produce  anomalous  subsurface  reflectors. 

b.  Increasing  impedance  environment.  The  initial  approach  assumed  a 
geologic  environment  of  increasing  impedance  layers.  In  other  words, 
sediments  become  more  competent  with  depth.  Whereas  this  assump¬ 
tion  holds  true  for  many  geologic  situations,  there  exist  numerous  situa¬ 
tions  where  soft  sediments  are  overlaid  by  more  competent  materials. 
The  AC50  program  presendy  does  not  incorporate  a  correction  for  this 
condition;  however,  as  discussed  in  detail  in  the  “Layer  Identification 
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Methods”  section  of  Qiapter  2,  techniques  are  provided  to  externally 
compensate  for  this  condition.  Do  not  rely  solely  on  the  AC50  pro¬ 
gram  for  sediment  characterization.  A  trial  version  of  the  AC50  pro¬ 
gram  that  performs  the  analysis  necessary  to  determine  the  direction  of 
the  impedance  change  is  now  being  fully  tested.  Initial  results  show 
significant  improvement  in  the  method’s  prediction  capabilities. 

c.  Natural  marine  sediments.  The  initial  impedance  function  used  for  this 
technique  is  based  on  empirical  data  collected  in  naturally  occurring 
marine  sediments  from  primarily  deeper  offshore  environments.  There¬ 
fore,  this  algorithm,  without  confirmatory  core  information,  may  pro¬ 
duce  anomalous  estimates  in  the  dynamic  nearshore,  harbor,  and 
riverine  environments-  the  primary  reason  for  the  regional  calibration 
approach. 

Whereas  these  assumptions  provide  a  practical  engineered  approach,  they 
may  also  limit  the  technique’s  ability  to  correctly  assess  a  specific  situation. 
There  is  no  one  single  sound  source  or  analysis  methodology  that  addresses  all 
engineering  and  geological  requirements. 


Limitations 

As  with  any  remote  sensing  technique,  limitations  exist.  These  limitations 
must  be  understood  in  order  to  appropriately  use  this  method.  Probably  the 
most  common  fault  encountered  in  geophysical  studies  is  improper  application 
of  a  given  technique  for  a  given  study  objective.  The  following  limitations 
exist  for  the  present  Acoustic  Impedance  technique. 

Signal-to-noise  ratio.  The  ability  to  accurately  assess  any  environment  is 
strictly  a  function  of  the  quality  of  data  obtained.  Low  S/N  data  will  pro¬ 
duce  poor  quality  results  or  possibly  no  results  at  all.  The  AI  method  limits 
its  processing  to  data  with  an  S/N  ratio  greater  than  5  db.  One  must  always 
be  suspicious  of  impedance  predictions  in  areas  of  poor  S/N.  Fortunately, 
most  noise  problems  can  be  corrected  through  effective  vessel  mobilization 
and  acoustic  calibrations. 

Layer  identification.  Unique  sediment  units  can  be  identified  only  when  an 
impedance  change  exists.  Gradual  changes  in  soil  type  may  not  result  in  an 
impedance  differential  large  enough  to  produce  a  reflection. 

Resolution.  Vertical  resolution  and  ultimate  depth  of  penetration  are  depen¬ 
dent  primarily  on  the  frequency  of  the  sound  wave.  Higher  operating  fre¬ 
quencies  permit  greater  resolution  of  marine  sediments  but  shallower  depths 
of  energy  penetration,  depending  on  the  characteristics  of  subbottom  mate¬ 
rials.  Also,  in  high-attenuation  sediments,  higher  frequencies  are  attenuated 
at  a  higher  rate  than  low  frequencies,  resulting  in  degradation  of  resolution 
and  errors  in  absorption  estimates  for  very  deep  layers.  As  previously 
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discussed,  a  number  of  field  techniques  and  processing  methodologies 
available  now  concurrently  improve  resolution  and  depth  of  penetration. 
However,  it  cannot  be  overemphasized  that  these  techniques  require  data 
with  high  S/N  ratios. 

Multiple  r^ections.  Multiples  are  one  of  the  primary  causes  of  data  qual¬ 
ity  degradation  in  shallow  marine  seismic  measurements.  Unlike  noise, 
which  can  be  distinguished  from  data  by  its  lack  of  lateral  continuity,  multi¬ 
ples  can  easily  be  mistaken  for  real  data,  create  false  structures,  and  change 
reflectivity  estimates.  Presently,  no  multiple-suppression  techniques  have 
been  developed  and  adapted  to  the  AI  system  of  software  products.  There¬ 
fore,  for  quantitative  an^ysis,  the  maximum  depth  of  investigation  is 
bounded  by  the  first  multiple  reflection  which  is  approximately  equal  to  the 
water  depth.  This  can  become  a  very  limiting  factor  in  shallow-water 
surveys. 

Beam  pattern  or  directivity.  Experience  has  shown  that  beam  pattern  and 
transducer  directivity  contribute  significanfiy  to  signal  degradation.  Sloping 
bottoms  and  rapidly  dipping  reflection  horizons  cause  inconsistent  reflection 
data  through  focusing  and  defocusing  of  the  incident  energy.  Rough,  irreg¬ 
ular  bottoms  with  numerous  scatterers  will  specularly  disperse  energy  away 
from  the  receiving  array.  For  quantitative  analysis,  minimization  of  the 
beam  angle,  either  through  beam  steering  or  receiver  focusing,  wlU  improve 
results  significantly  and  will  likely  be  incorporated  in  the  WES  AI  system. 
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6  Conclusions  and 
Recommendations 


Conclusions 

In  its  present  state  of  development,  acoustic  impedance  processing  of  seis¬ 
mic  reflection  data  provides  a  reasonably  accurate,  continuous  description  of 
bottom  and  subbottom  marine  sediment  characteristics  in  a  rapid,  cost-effective 
manner.  From  a  statistical  standpoint,  density  can  be  acoustically  derived  to 
within  ±10  percent  of  in  situ  conditions.  Properly  calibrated  surveys  provided 
Corps  Districts  and  dredging  contractors  with  the  following  results: 

a.  Density  estimates  of  marine  sediments. 

b.  Continuous  subbottom  information  for  planning  and  designing  dredging 
and  sampling  programs. 

c.  Estimates  of  the  volume  and  type  of  material  to  be  removed  through 
dredging. 

d.  A  detailed  and  continuous  geologic  database  for  aiding  long-tenn  plan¬ 
ning  of  future  work. 

Acoustic  impedance  information,  if  properly  implemented  in  the  project 
planning  stages,  provides  valuable  data  on  the  distribution  and  extent  of  dif¬ 
fering  marine  sediments,  aids  in  locating  optimal  placements  of  sampling 
cores,  and  supplements  previously  obtained  soil  borings  by  providing  continu¬ 
ous  profile  coverage  of  sediment  characteristics  between  sample  locations. 

Further,  the  acoustic  impedance  method  is  not  limited  to  sediment  charac¬ 
terization  but  can  also  provide  such  essential  infonnation  as: 

a.  Location  of  marine  sand  deposits  for  beach  replenishment. 

b.  Long-term  monitoring  of  dredged  material  disposal  areas. 

c.  Delineation  of  submarine  geologic  fonnations. 
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d.  Detection  of  submarine  features  such  as  pipelines  or  other  dredging 
hazards. 

e.  Identification  of  fluid  mud  in  navigable  waterways. 


Recommendations 


It  is  recommended  that  development  continue  beyond  the  DRP  in  expand¬ 
ing  the  capabilities  of  the  waterborne  geophysical  exploration  system.  Empha¬ 
sis  should  be  placed  on  expansion  of  capabilities  to  include  a  wide  range  of 
geophysical  techniques.  Specifics  include  the  following: 


Hardware 

a.  Multichannel  seismic  acquisition  and  processing  system.  The  benefits 
of  multichaimel  processing  include  S/N  improvement,  multiple  and  dif¬ 
fraction  suppression,  increased  knowledge  of  the  subbottom  velocity 
profile,  and  the  potential  use  of  2-D  algorithms. 

b.  Automated  gain  recording.  The  acquisition  system  should  be  upgraded 
to  digitally  record  any  amplifier  gain  changes  occurring  during  data 
acquisition.  This  relatively  simple  upgrade  will  eliminate  potential 
errors  in  recording  field  data  and,  more  importantly,  will  allow  for 
maximum  flexibility  in  acquiring  highest  S/N  data. 

c.  Minimize  beam  angle.  For  quantitative  analysis,  minimization  of  the 
beam  angle  will  improve  results  significantly.  This  can  be  accom¬ 
plished  either  at  the  source  or  receiver.  Source  beam  angle  is  a  func¬ 
tion  of  transducer  design.  Narrow-beam  transducers  for  the  low 
operating  frequencies  (<10  deg)  are  generally  very  large  arrays,  in  most 
cases  impractical  for  small-boat,  shallow-water  applications.  Specially 
designed  receiving  hydrophone  arrays  can  be  constructed  and  used  in 
conjunction  with  “beam-forming”  software  to  more  practically  accom¬ 
plish  the  task.  Such  systems  are  commercially  available  and  could  be 
readily  integrated  with  the  existing  system. 

d.  Core  analysis  standardization.  Minimum  core  analysis  requirements 
should  be  implemented  to  specifically  address  geoacoustic  analysis. 
Minimum  requirements  should  be  as  follows: 

(1)  Independent  direct  measurement  of  acoustic  attenuation  of  the 
sediment  material  retrieved  from  the  seafloor. 

(2)  Direct  measurement  of  compressional  (p-wave)  velocity  of  core 
sediment. 

(3)  Determination  of  specific  gravity,  density,  and  porosity. 
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(4)  Grain  size  analysis  of  all  sediment  types,  from  fine-  to  coarse¬ 
grained  sediments. 


Software 

e.  Dependable  algorithm  for  determining  sign  of  reflection  coefficient. 

The  upgrade  is  highly  recommended  for  the  purpose  of  assessing 
decreasing  impedances  in  successive  layers.  This  algorithm  should  be 
applied  directly  in  the  AC50  program.  An  initial  AC60  program  is 
now  undergoing  beta  testing. 

/.  Multiple  suppression  algorithm.  This  algorithm  is  needed  to  suppress 
false  reflectors  and  provide  better  clarification  of  real  reflectors,  leading 
to  better  estimates  of  reflectivity,  impedance,  and  other  geotechnical 
parameters.  It  would  also  increase  the  effective  depth  of  penetration  in 
shallow-water  studies  as  well  as  improve  the  reliability  of  deeper 
reflectors. 

g.  Improved  forward  modeling  algorithms.  More  comprehensive  methods 
for  generating  synthetic  seismograms  of  geologic  data  are  needed  to 
better  model  absorption  and  internal  reflections. 


Research  needs 

h.  Geoacoustic  relationships.  An  extensive  database  of  acoustic  and  geo¬ 
technical  data  has  been  acquired  during  the  development  of  this 
method.  Detailed  study  into  the  geoacoustic  relationships  of  nearshore 
sediments  is  needed,  specifically: 

(1)  Grain  size  analysis  (sorting,  skewness,  kurtosis,  etc.)  for  the  pur¬ 
pose  of  better  describing  sediments  for  applications  such  as  beach 
nourishment  projects. 

(2)  Analysis  of  contaminated  sediments. 

(3)  Analysis  of  highly  fluidized  sediments  such  as  fluff/fluid  mud. 

i.  Soil  properties.  The  database  should  be  utilized  to  study  the  relation¬ 
ships  between  different  soil  properties  such  as  density/porosity  versus 
mean  grain  size.  This  will  aid  in  developing  better  sediment  charac¬ 
terizations  using  geoacoustic  methods. 
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Rgure  68.  Photographs  of  WATERWAYS  EXPLORER  (Sheet  1  of  3) 


Figure  68.  (Sheet  2  of  3) 
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Figure  68.  (Sheet  3  of  3) 
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Appendix  A 
Equipment 


Hardware 


The  study  of  the  AI  method  did  not  attempt  to  develop  and  introduce  new 
acoustic  sound  sources  or  receivers.  Rather,  it  was  decided  to  utilize  the 
extensive  array  of  commercially  available  systems,  making  the  Al  method  both 
flexible  and  robust  enough  to  meet  the  challenge  of  assessing  a  wide  variety  of 
marine  environments  for  commercial  and  government  surveyors. 

A  wide  range  of  shallow  seismic  data  acquisition  systems,  beginning  with 
analog  and  transitioning  to  digital,  were  utilized  during  the  development  of  the 
AI  method.  It  is  important  to  emphasize  the  fact  that  a  variety  of  off-the-shelf 
commercially  available  equipment  can  be  adapted  to  AI  technology.  A 
detaUed  description  of  every  system  is  not  appropriate  within  this  text.  A 
block  diagram  showing  equipment  used  in  the  WES  integrated  geophysical 
system  is  shown  in  Figure  Al.  This  system  is  currently  used  to  perform  AI 
surveys  and  other  types  of  marine  geophysical  investigations.  Specific  systems 
utilized  for  the  AI  application  are  as  follows: 

a.  J.5-  7.0-kHz  pinger  system.  This  system  is  by  far  the  most  f)opular 
utilized  in  the  shallow  marine  seismic  industry.  It  allows  transmission 
of  variable-length  pulses  (0.2  -  3  msec)  of  3.5-  and  7.0-kHz  frequen¬ 
cies.  Power  levels  can  be  varied  from  1  to  10  kW.  Figure  A2  shows 
an  over-the-side  deployment  of  a  pinger  system.  However,  depth  of 
penetration  can  be  limited  in  areas  of  highly  competent  (dense)  sedi¬ 
ments.  Considerable  efforts  were  expended  in  developing  alternative 
deployment  configurations  to  improve  S/N  characteristics  of  pinger  data 
in  the  shallow  marine  environment  These  systems  were  originally 
designed  to  operate  in  water  depths  greater  than  about  50  ft,  resulting  in 
configurations  employing  integrated  transmit/receive  (T/R)  networks  in 
order  to  utilize  the  same  transducers  as  transmitters  and  receivers.  The 
resulting  transducer  ringing  and  coupling  creates  coherent  noise  keyed 
to  the  transmitter  timing.  In  shallow  water  (less  than  30  ft),  significant 
S/N  problems  arise  due  to  coherent  noise  from  the  transmitter  inter¬ 
fering  with  the  first  return.  Figure  A3  illustrates  this  problem. 
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To  improve  the  S/N,  a  separate  receiving  array  is  now  deployed  indepen¬ 
dent  of  the  transmitter.  By  decoupling  the  receiving  array  from  the  transmitter 
and  physically  separating  the  transducers,  all  of  the  near-field  transmitter  ring¬ 
ing  is  eliminated  from  the  bottom  reflection,  regardless  of  water  depth.  Fig¬ 
ure  A4  presents  reflection  data  collected  with  this  configuration.  This  has 
become  the  standard  pinger  and  chirper  deployment  configuration  for  the  AI 
method. 

b.  Boomer  system.  This  system  is  a  high-energy,  medium-bandwidth  unit 
providing  up  to  1,000  J  of  energy  in  the  400-  to  5,000-Hz  frequency 
band.  The  system  is  designed  to  provide  reasonable  vertical  resolution 
combined  wifli  greater  penetration  depths  in  more  competent  sediments. 
Because  of  the  high  power  involved  and  because  the  coherent  noise 
radiates  to  the  receiver  as  well  as  to  the  bottom,  a  separate  towed 
receiving  array  is  used.  This  array  is  normally  towed  at  right  angles  or 
directly  aft  of  the  source.  The  exact  tow  point  is  detennined  by  the 
water  depth  and  by  minimizing  the  coherent  noise  to  the  receiver.  In 
practice,  the  source-to-receiver  separation  can  be  fine-tuned  by  an 
experienced  data  observer  to  produce  the  best  record  quality  obtainable 
at  a  specific  site.  Figure  A5  shows  a  typical  boomer  deployment. 

c.  Chirp  system.  Commercially  available  chirp  systems  (analogous  to  land 
surface  vibrators  used  in  the  petroleum  industry)  are  designed  to 
improve  the  S/N  and  vertical  resolution  by  the  application  of  correlation 
processing  (matched  filtering)  to  a  wide-band  swept-frequency  long- 
pulse  signal.  The  improvement  is  a  function  of  the  bandwidth-time 
product  Standard  chirp  systems  normally  operate  from  1  kHz  to 

10  kHz  with  pulse  lengths  of  3  to  20  msec.  At  this  point,  it  must  be 
stated  that  only  the  chirper  transducer,  not  a  complete  “off-the-shelf’ 
chirper  data  acquisition  system,  was  used  in  development  of  the 
AI  method.  Correlation  processing  of  chirp  data  will  be  discussed  in 
detail  later  in  the  software  section  of  this  chapter.  Since  the  mid¬ 
frequency  energy  range  of  most  chirp  systems  is  around  4  kHz,  the 
depth  of  penetration  achieved  is  about  the  same  as  pinger  systems.  As 
with  the  pinger  system,  the  receiving  array  is  deployed  independently  of 
the  transmitting  array.  This  procedure  greatly  improves  the  S/N  by 
eliminating  the  coherent  noise  resulting  from  transmitter  ringing.  The 
complete  chirp  data  acquisition  system  developed  by  WES  allows  for 
unprecedented  flexibility  in  generating  source  wave  fonns.  In  many 
cases,  the  need  arises  for  tailoring  the  seismic  source  to  produce  maxi¬ 
mum  energy  within  a  predetermined  frequency/wavelength  band  for 
resolution  of  a  site-specific  target.  Chirper  transducer  control  software 
was  developed  by  WES  to  allow  the  design  of  virtually  any  wave  form 
type  (within  the  design  limits  of  transducer  response)  as  a  function  of 
frequency,  amplitude,  and  pulse  length.  This  wave  fonn  is  delivered 
directly  to  the  power  amplifier  of  tlie  chirp  transceiver  through  a  digi- 
tal-to-analog  signal  converter,  linear  amplifier,  and  analog  filter. 


A2 
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d.  Bubble  pulse  system.  The  bubble  pulser  generates  a  low-  to  mid-range 
frequency  wavelet,  with  a  frequency  content  between  400  and 
2,000  Hz,  with  most  of  the  energy  concentrated  between  600  and 
900  Hz.  Because  of  the  source’s  low-frequency  content,  penetration 
depth  in  competent  materials  such  as  sands  is  signiflcantly  greater  than 
with  the  3,5-kHz  system.  However,  increased  depth  of  penetration 
comes  at  the  cost  of  resolution.  Because  the  bubble  pulser  exhibits 
bandwidth  rather  than  single-frequency  energy,  the  same  correlation 
techniques  to  improve  resolution  in  chiip  technology  can  be  applied, 
improving  the  and  resolution  by  a  factor  of  two. 

Based  upon  the  results  of  numerous  surveys,  it  is  recommended  that  multi¬ 
ple  systems,  preferably  with  different  frequencies  and  energy  levels,  be  utilized 
for  all  AI  surveys.  For  most  Corps  dredging  applications,  high-resolution 
delineation  of  surface  sediment  layers  is  important,  requiring  higher  frequency 
devices  and/or  chirp  technology.  In  order  to  accurately  assess  the  absorption 
characteristics  of  subbottom  sediments  (required  to  characterize  the  subbottom 
sediments  discussed  in  Chapter  2,  ‘Technical  Development”)  multiple  frequen¬ 
cies  encompassing  the  greatest  bandwidth  are  required. 


Software 

Reflected  acoustic  signals  detected  by  the  receiving  array  are  first  amplified 
and  filtered  as  needed  to  ensure  maximum  S/N.  The  amplifiers  must  be  linear 
and  must  exhibit  no  DC  bias.  The  AI  method  requires  precise  knowledge  of 
the  total  energy  due  to  signal  amplification.  Real-time  filtering  is  provided  at 
the  front  end  of  amplifiers  to  reduce  undesirable  noise  as  much  as  possible 
prior  to  digital  sampling  of  data.  After  this  preprocessing,  data  are  recorded 
digitally  using  a  specially  designed  Digital  Field  Shallow  Seismic  Acquisition 
System  (DFAS).  A  complete  description  and  detailed  discussion  of  the  DFAS 
operation  is  provided  by  Caulfield  (1991a). 

The  DFAS  is  designed  to  provide  an  economical  means  of  recording  shal¬ 
low  marine  seismic  data  on  commonly  available  computing  systems.  The 
DFAS  is  an  IBM-compatible  hardware/software  package  that  operates  under 
DOS  3.3  or  greater.  The  system  has  a  minimum  dynamic  range  of  72  db  and 
provides  real-time  visual  color  display,  disk  writing  procedures,  and  a  data 
processing  and  playback  system. 

Real-time  data  acquisition  is  based  on  the  data  input  being  wide  band  and 
having  sufficient  S/N  levels  to  meet  standard  communications  quality  control 
conditions.  Acquisition  quality  control  is  provided  to  verify  that  the  amplifier 
gains  are  set  for  optimum  S/N,  appropriate  sample  rates  and  trace  lengths  are 
chosen,  and  appropriate  timing  offsets  are  employed  to  maximize  the  reflection 
window.  The  system  is  designed  to  handle  most  standard  shallow  subbottom 
geophysical  tools,  such  as  3.5-kHz  pingers,  boomers,  etc.  The  analog/digital 
converter  has  a  12-bit,  20-msec  sampling  rate  with  precision  sample  and  hold 


Appendix  A  Equipment 


A3 


amplifiers  on  the  front  end.  The  actual  operational  sample  rate  is  dependent 
on  the  host  computer  and  is  software  controlled. 

The  system  has  been  designed  to  operate  with  any  user-supplied  navigation 
system  that  has  an  RS-232  output.  Navigation  information  is  read  directly  into 
data  file  headers  for  direct  correlation  with  subbottom  data. 


Limited  postprocessing  options  are  provided.  These  include  horizontal 
spatial  stacking,  noise  reduction,  and  linear  gain  manipulation.  Also,  data  may 
be  played  back  utilizing  the  spherical  spreading  correction  to  compensate  for 
transmission  loss.  Typical  data  records  for  pinger,  boomer,  and  bubble  pulser 
systems  are  presented  in  Figures  A6,  A7,  and  A8,  respectively. 

Chirp  data  are  recorded  using  Real-Time  Correlation  Acquisition  (RTCIO) 
System  shallow  seismic  software.  This  system  provides  real-time  matched 
filter  correlation  processing  of  echo  time  series  with  a  true  replica  of  the  out¬ 
going  source  wavelet  The  match  filter  correlation  process  utilized  here  is 
described  in  Chapter  2  (‘Technical  Development”),  and  a  detailed  description 
of  the  RTCIO  program  is  provided  by  Caulfield  (1993). 


Figure  A1.  Block  diagram  of  WES  integrated  geophysical  system 
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iure  A3.  Coherent  noise  due  to  couplirvg  of  transmitter  and  receiver 
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Figure  A4.  Improvement  in  signal-to-noise  due  to  separated  configuration 


Figure  A5.  Typical  boomer  system  in  tow 
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Figure  A6.  Typical  pinger  data  record 


Figure  A7.  Typical  boomer  data  record 


Number 
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Figure  A8.  Typical  bubble  pulse  data  record 
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US  Army  Corps 
of  Engineers 

New  England  Divisbn 


Boston  Harbor,  MA 
Navigation  improvement  Project 
Information  Sheet 

Uuly  1992 


Feasibility  Report  -  Final  Feasibility  Report  and  Environmental  Assessment  completed 
in  May  1988  in  response  to  Senate  resolutions.  Committee  on  Public  Works,  1  March  1968 
and  11  September  1969.  Report  of  Board  of  Engineers  for  Rivers  and  Harbors  (BERH)  was 
completed  on  21  September  1988.  The  Chief  of  Engineers  Final  Report  is  dated  11  May 
1989.  The  project  was  authorized  in  the  Water  Resources  Development  Act  of  1990,  signed 
on  Novemter  28,  1990. 

While  Boston's  principal  entrance  and  main  access  channels  are  40  feet  deep,  the  3  major 
tributaries,  along  which  the  majority  of  port  terminals  are  located,  are  35  feet  deep.  The 
resulting  tidal  delays  and  limits  on  vessel  size  and  loading  place  significant  restrictions  on 
shippers.  The  objective  of  the  study  was  to  substantially  reduce  the  cost  of  transporting 
bulk  commodities  through  the  Port  of  Boston  by  eliminating  or  reducing  tidal  delays  for 
larger  vessels. 

Alternative  plans  investigated  included  deepening  of  each  of  the  three  tributaries,  the 
lower  Mystic  River  (currently  authorized  as  a  separate  project),  Chelsea  River  and 
Reserved  Channel  were  examined  at  incremental  depths  up  to  42  feet.  Also  a 
nonstructural  improvement  entailing  designation  of  specific  Federal  channel  limits 
through  the  middle  reaches  of  the  harbor  would  result  in  an  enlarged  deepwater 
anchorage,  at  no  cost. 


Project  Description  (See  Attached  Figure): 

•  Reserved  Channel  -  Narrow  the  35-foot  channel  from  430  to  400  feet  wide,  widen  the 
channel  at  its  confluence  with  the  Main  Ship  Channel,  relocate  Harbor  lines  at  the 
confluence  and  deepen  the  channel  to  40  feet,  except  for  its  upper  1^40  feet  which  'would 
remain  at  35  feet.  Lx>cal  dredging  of  widened  berthing  areas  along  the  southern  limit  to  40 
feet  or  greater. 

•  Chelsea  River  Channel  -  Deepening  the  existing  35-foot  channel  to  38  feet.  Local 
dredging  of  berthing  areas  at  the  major  petroleum  terminals  to  38  feet  or  greater  and  local 
relocation  and  alteration  of  .utility  crossings  beneath  the  channel. 


Figure  B1.  Boston  Harbor,  MA,  Navigation  Improvement  Project  information  sheet 
(Sheet  1  of  6) 
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•  Mystic  River  Project  -  Deepen  the  existing  35-foot  channel  to  40  feet,  except  for  an  area  at 
the  upstream  limit  along  the  southern  shoreline  where  the  waterfront  has  ^n  converted 
to  non-navigation  dependent  uses  and  where  existing  lisers  do  not  require  depths  greater 
than  35  feet.  Local  dredging  of  major  terminal  berthing  areas  to  40  feet  or  greater. 

•  Tnnpr  Ponfluence  Area  -  Widen  the  35-foot  confluence  of  the  Mystic  and  Chelsea  Rivers 
and  deepen  to  40  feet  to  enable  deepening  of  the  two  river  channels. 

•  President  Roads  Ship  Channel  -  Nonstructural  designation  and  remarking  of  specific 
Federal  channel  limits  along  the  southern  reach  of  the  Roads  connecting  the 
confluence  of  the  three  entrance  channels  which  the  Inner  Harbor  Main  Ship  Channel, 
resulting  in  a  20%  enlargement  of  the  deepwater  anchorage  to  about  420  acres. 

Pertinent  Data  on  Project  Features: 

.  Channel  Dredging  -  2.038.000  cy  +  94.000  cy  rock  =  2.132.000  cy 

•  Berthing  Area  Dredging  by  non-Federal  interests  (to  -  40  feet)  -  228,700  cy 

•  Utility  Relocations  -  Chelsea  River  Channel  only  -  tower  Boston  Edison  aixJ  MBTA  electnc 

cables,  arxl  remove  abandoned  MWRA  water  tunnel.  Provide  increased 
protection  to  Boston  Gas  gas  siphon. 

•  Channel  Limit  Designation/Realignment  -  President  Roads  Ship  Channel  -  Designate 

specific  Federal  channel  limits  through  the  Roads  to  jurisdictionally  connect 
entrance  channels  with  the  Main  Ship  Channel  and  enlarge  the  anchorage  area 
by  70  acres  in  deepwater. 

views  of  Others;  The  Massachusetts  Port  Authority  (Massport),  the  sponsor,  has 
indicated  its  willingness  and  capability  to  act  as  project  sponsor  stated  iii  a  14  March 
1989  letter.  Harbor  users,  local  chanters  of  commerce,  etc.,  are  supportive  of  the 
prpject. 

The  Secretary  of  Environmental  Affairs  provided  a  letter  of  strong  support  for  this 
project  dated  14  December  1989.  State  agencies  condition  their  supjwrt  for  continued 
study  on  the  suitability  of  dredging  and  ocean  disposal  of  the  material. 

Federal  agendes  condition  their  approval  of  the  project  on  the  design  of  an  acceptable 
dredge  disposal  plan. 

The  U.S.  Coast  Guard  expressed  concern  for  navigational  safety  in  the  Chelsea  River 
Channel  in  the  vidnity  of  the  Chelsea  Street  Bridge.  The  City  of  Boston  plans  to  replace 
the  dilapidated  fender  system. 
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Cost  of  the  Recommended  Project: 

Estimated  Implementation  Costs:  Fully  Funded  Estimate  (includes  inflation) 

Federal  (Commercial  Nav^tion) 

Corps  -  65%  of  General  Nav.  Features  (GNF)  Cost 
U.S.  Coast  Guard  -  Rebcation  Nav.  Aids 
TOTAL 

Non-Federal 

25%  Of  GNF  Cost  -  Up  Front 
1 0%  GNF  Cost  -  Payback  (minus  reloc.  credX) 

1 00%  Berthing  Area  Deepening 
Utility  Relocation  -  Chelsea  River  Oniy 
TOTAL 

TOTAL  PROJECT  COST 

Project  Economic  Feasibility  (March  1992); 

Project  Economic  Life  -  50  years,  Base  Year  1995 
Benefit  -  Cost  Ratio:  1.68  (Discount  Rate  FY92  8-1/2) 

Project  Effects: 

Environmental  Quality;  Consfruction  of  the  project  will  be  scheduled  in  such  a  way  as 
to  avoid  adverse  impacts  on  aquatic  life.  Analysis  of  environmental  testing  of  dredge 
material  has  found  that  material  from  the  Reserved  channel  is  suitable  for  open  ocean 
disposal  and  that  sediments  from  the  Mystic  and  Chelsea  channels  are  unsuitable  for 
unconfined  ocean  disposal.  Capping,  using  suitable  project  material,  is  included  in  the 
cost  estimate  and  schedule. 

Other  Effect^.  As  the  New  England  Region's  largest  and  most  active  port,  shipping  over 
20  million  tons  of  cargo  annually,  improvements  to  Boston  Harbor  will  enhance  the 
regional  economy  through  lowering  transportation  costs  for  bulk  and  containerized 
cargo.  The  deepening  of  existing  channels  using  a  waterborne  plant  will  not  result  in 
adverse  social  impacts. 


$21,200,000 

35.000 

$21,215,000 


$  7,790.000 
2.160,000 
2,400.000 
950.000 
$  13,300,000 
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The  Reserved  Channel: 

Coastal  Oil  New  England  Inc.(Petroleum),  Boston  Edison  (Petroleum  for  Power 
Generation),  and  shipping  lines  using  Massport’s  Conley  Container  Terminal. 


The  Mystic  River  Channel: 

Exxon  USA,  lnc.(Pefroleum),  Distrigas  of  New  England  (LNG) 

Prolerized  (Scrap  Metal  Export),  Boston  Edison  (Petroleum  for  Power  Generatiohj 

All  beneficiaries  on  the  Chelsea  River  Channel  are  petroleum: 

Moba  Oil,  Global ,  Northeast  Petroleum,  Gulf  OO  Co.,  and  Coastal  Oil  New 
England  Inc.. 

Harbor  users  along  the  project  channels  have  contacted  to  reevaluate  project  benefits. 
Users,  not  currently  listed  above  as  direct  beneficiaries  will  be  added  if  it  can  be  shown 
that  the  project  will  improve  their  operations.  The  results  of  this  analysis  should  be 
available  in  July  1992 

Project  Implementation  Schedule: 


Initiate  Design 

September  1990 

Design  Memorandum  App^rv’d 

May  1993 

Sign  Local  Cooperation  Agreement 

November  1994 

Complete  Plans  and  Specificaticns 

February  1994 

Award  Construction  Contract 

February  1995 

Project  Completion 

May  1996 
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other  Information: 

Sponsor  Supporting  Activity  -  Massport  is  actively  working  with  terminal  owners  for 
joining  forces  to  perform  environmental  testing  for  dredge  material  to  be  removed 
from  berthing  areas. 

Maissport  will  obtain  the  required  permits  from  Federal  and  state  agencies.  They  began 
this  process  on  8  May  1991  by  filing  an  Environmental  Notification  Form  (ENF)  with 
the  MEPA  Office.  The  Secretary  of  Environmental  Affairs  Certificate  was  signed  7  June 
1991,  An  Environmental  Impact  Report  (EIR)  is  under  preparation  by  Massport  through 
a  contract  with  Normandeau  Associates.  The  final  permits  must  be  obtained  prior  to 
the  signing  of  the  Local  Cooperation  Agreement  sdieduled  for  November  1994. 

EIS  -  The  Corps  decided,  in  April  1992,  to  prepare  an  EIS  for  the  project.  This  decision 
was  based  on  the  substantial  diange  in  the  volume  of  material  to  be  dredged  (561,000  cy 
of  material  from  berths  that  are  not  project  beneficiaries),  and  the  significant  public  and 
agency  controversy  regarding  the  capping  alternative  for  disposal.  ITie  EIS  process  will 
parallel  that  of  the  Massport  EIR  and  will  result  in  a  combined  EIS/EER  document. 

There  is  no  direct  impact  on  the  environmental  impact  phase  of  the  project  design, 
however,  the  process  will  delay  the  request  for  a  new  construction  start  from  June  1992 
to  June  1993.  The  impact  of  this  decision  are  included  in  the  costs  and  schedule  shown 
herein. 

Project  Management  Plan  -  The  Project  Management  Plan  (PMP)  has  been  approved  by 
NED  and  Massport.  Modifications,  to  include  changes  made  since  August  1991  and 
summarized  herein,  are  being  prepared. 

Chelsea  River  Issues  -  Massport  sponsored  two  workshops  for  Chelsea  Creek  interests 
to-discuss  issues  related  to  petroleum  terminals.  These  "tank  farms"  supply  one  quarter 
of  the  regional  heating  oil  and  over  two  thirds  of  the  gasoline  as  well  as  ofiier  products. 
The  narrow  opening  through  the  Chelsea  Street  Bridge  severely  limits  the  size  of 
tankers  using  upstream  terminals.  The  collapse  of  a  bulkhead  wall  and  blockage  of  the 
channel  in  January  1990  add  to  the  vulnerability  of  the  region  to  Chelsea  River 
shipping,  storage  and  distribution.  NED  has  participated  in  these  workshops,  "nie 
recommended  project  will  improve  present  conditions  but  can  not  resolve  all  issues 
facing  this  area.  There  has  also  been  renewed  interest  in  the  replacement  of  the  Chelsea 
Street  Bridge  by  the  Qty  of  Boston  and  the  Coast  Guard.  The  Coast  Guard  is  proceeding 
with  the  required  efforts  to  obtain  funds  for  a  new  bridge  through  the  Truman-Hobbs 
Act. 
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Ship  Simulation  Model  Studies  -  A  hydrodynamic  model,  to  determine  changes  in 
currents  caused  by  the  project,  has  been  developed  at  NED.  Results  from  this  model 
has  been  used  in  a  ship  simulation  model  whiA  is  required  to  assure  that  the  project  is 
both  safe  and  efficient  A  deeper  channel,  which  changes  currents,  and  heavier  ships 
will  effect  the  response  and  handling  of  the  pilots.  Ship  Simulation  is  underway  at 
Marine  Safety  International  (MSI),  in  Newport  RI  and  Kings  Point  NY.  MSI  initiated 
work  on  14  October  1991  and  is  scheduled  for  completion  in  October  1992.  Verification 
runs  of  existing  channel  conditions  were  completed  in  February  1992.  Tests  for  project 
conditions  were  completed  in  April  1992.  Analysis  of  the  results  which  could  lead  to 
design  changes,  if  any,  will  be  completed  in  July  199Z 

Disposal  of  Dredged  Material  -  Sediments  dredged  from  the  Mystic  and  Chelsea 
channels  contain  some  contamination  which  makes  it  unsuitable  for  unrestricted 
ocean  disposal.  The  project  cost  and  schedule  includes  capping  the  contaminated 
sediments  with  day  which  imderlies  the  contaminated  material.  There  is  suffident 
volume  for  the  day  in  each  channel  to  provide  an  adequate  cap  for  it's  sediments. 
Because  the  MBDS  is  located  in  deep  water  (80  to  100  meters),  there  is  concern  from 
both  Federal  and  state  agendes  that  capping  may  not  be  effective  in  protecting  the  water 
from  contamination.  The  Corps  has  done  extensive  research  and  pilot  studies  to 
develop  techniques  for  capping.  A  workshop  held  on  15  October  1991  brought  together 
all  concerned  agendes  to  discuss  disposal  issues.  Their  concerns  will  be  addressed  in  a 
dredge  material  disposal  plan  which  will  indude  design,  construction  and  monitoring  . 

Disposal  at  MBDS  may  also  be  subject  to  more  extensive  controls  depending  on  the 
final  plan  for  the  proposed  Stellwagen  Bank  Sanctuary. 

Economic  Re-evaluation  -  During  July  1992  results  of  an  economic  survey  of  project 
benefidaries  will  be  completed.  The  objective  of  this  survey,  and  the  analysis  which 
will  follow,  is  to  update  the  economic  justification  of  the  project  last  done  in  1988. 
Benefidaries  may  be  added  to  the  project  if  their  operations  become  more  efficient  from 
the  deeper  channels. 

Other  Major  Projects  and  Issues  -  There  are  two  major  projects  underway  in  Boston 
Harbor  which  will  continue  during  the  same  timeframe  as  the  Navigation 
Improvement  Project:  MWRA  projects  for  the  deanup  of  the  harbor  and  the  Central 
Artery  /  Third  Itobor  Crossing.  Close  coordination  with  these  projects  will  be  essential. 

There  is  major  concern  for  the  extent  and  toxidty  of  the  industrial  waste  site  adjacent  to 
the  MBDS.  The  EPA  is  investigating  this  issue  and  may  recommend  solutions.  Close 
coordination  with  EPA's  work  will  continue. 
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CEWES-GG  15  April  1992 

MEMO  FOR:  Thomas  R.  Patin  (EP-D) 

FROM:  ROBERT  F.  BALLARD,  JR. 

SUBJECT:  DOTS-Sponsored  Trip  to  New  England  Division  (NED) 


1.  As  a  result  of  a  letter  request  from  NED  dated  6  February 
1992,  Subject:  Dredging  Operations  Technical  Support,  a  technolo¬ 
gy  transfer  visit  was  made  to  NED  on  7  April  1992.  Messers. 
Richard  G.  McGee  (HS-H) ,  Keith  J.  Sjostrom  (GG-F) ,  Thomas  S. 
Harmon  (GG-F) ,  and  Robert  F.  Ballard,  Jr.  (GG)  were  the  WES 
personnel  in  attendance  at  the  meeting.  The  purpose  of  the 
meeting  was  for  WES  personnel  to  present  a  seminar  on  the  appli¬ 
cation  of  Acoustic  Impedance  Technology  to  the  upcoming  Boston 
Harbor  Dredging  Improvement  Project. 

2.  Mr.  Peter  Jackson,  Project  Manager  of  the  Boston  Harbor 
Project  was  point  of  contact.  Mr.  Jackson’s  request  for  the 
seminar  came  about  as  the  result  of  his  having  seen  the  WES 
Dredging  Research  Program  video  "Acoustic  Impedance  Technology" 
produced  by  DRP  work  unit  #32470  "Rapid  Measurement  of  Properties 
of  Consolidated  Sediments". 

3.  Prior  to  the  meeting  at  NED,  Mr.  Jackson  sent  me  sufficient 
background  information  on  the  overall  project  to  enable  us  to 
plan  a  presentation  tailored  specifically  to  his  application. ^  It 
was  felt  that  four  WES  people  should  participate  in  this  meeting 
because  each  is  a  specialist  and  could  perform  as  an  expert 
during  a  formal  tutoral  presentation  and  answer  questions  posed 
in  open  forum  at  the  end  of  the  seminar.  Approximately  20  people 
were  in  attendance  during  the  meeting.  Representatives  in 
addition  to  NED  personnel,  were  from  the  Massachusetts  Port 
Authority,  USGS,  and  the  EPA. 


4.  The  format  and  content  of  the  seminar  can  be  gleaned  from 
the  attached  hard  copies  of  my  Vu-Graphs.  Our  reception  was  both 
hospitable  and  communicative.  Audience  participation  was  excep¬ 
tional.  We  were  told  that  an  intelligent  decision  could  now  be 
made  by  NED  regarding  the  possible  use  of  Acoustic  Impedance 
during  the  planning  phase  of  the  Boston  Harbor  Project.  We  will 
be  contacted  by  Mr.  Jackson  when  our  degree  of  participation  is 
decided  upon. 

5.  Closing  comment:  I  am  under  the  impression  that  ^his  under¬ 
taking  is  a  prime  example  of  excellent  cooperation  between  the 
complementary  missions  of  DOTS  and  the  DRP. 

cc:  Dr.  Don  Banks  (GS) 

Mr.  Clark  McNair  (CP-D) 

Mr.  Russ  Tillman  (CP-D) 
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To  Acquaint  You.  the  Potential  User, 
With  the  Acoustic  Impedance 
Procedure 

Development  Background 
Advantages  and  Limitations 
Application  to  Your  Project 
Benefits  vs.  Cost 
Spinoff  Applications 
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■  Spinoff  Possibilities 
-Anomaly  Detection 
-Utility  Locations 
-Bottom  Surveillance 


Project 

■Axxist»  Impedance  De«ifeiopment  Background 
— Mvantages  and  UmXatiocis 

■  Discussion  of  Application  to  NED  Projects 

-Six  Passible  Sites 

■  Recommendations 
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-Side  Scan  Sonar 

■  Support  Requiiemenls 

-People  and  Equipment 
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When  Properly  Applied  During 
the  Planning  Stages  of  a 
Major  Project,  the  Acoustic 
Impedance  Survey  Will 
Appreciably  Reduce  Costs  And 
Provide  a  Database  for 
Perpetuity 
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SUBBOTTOM  HYDRO-ACOUSTIC  SURVEY  OF  THE 
BOSTON  HARBOR  TRIBUTARIES:  RESERVED  CHANNEL; 
MYSTIC  RIVER;  CHELSEA  RIVER;  AND,  INNER  CONFLUENCE 


Background 

1.  A  design  memorandum  for  improvement  dredging  of  the  three 
tributaries  to  Boston  Harbor  are  soon  to  be  prepared  by  the  New 
England  Division.  A  description  of  the  materials  to  be  dredged 
in  terms  of  density  and  soil  classification  are  necessary  for 
planning  and  monitoring  of  the  dredging  program  and  disposal. 

The  USAE  Waterways  Experiment  Station  (WES)  is  involved,  through 
the  Dredging  Research  Program  (DRP) ,  in  developing  a  rapid 
geophysical  technique  to  determine  material  characteristics  of 
bottom  and  siabbottom  sediments.  Acoustic  impedance  values 
determined  from  the  seismic  reflection  data  are  directly  related 
to  the  density  and  material  type  of  the  subbottom  sediments 
which  in  turn  are  factors  related  to  the  degree  of  difficulty  in 
removing  this  material  through  dredging.  Results  to  date 
indicate  good  correlation  between  the  values  determined  with 
this  technique  and  insitu  information  collected  at  several 
sites. 

Objective 

2.  The  project  objective  is  to  quantify  with  depth  the  bottom 
and  subbottom  sediments  in  terms  of  material  density  and  soil 
classifications  to  a  minimum  depth  of  -42  feet  Mean  Low  Water 
(MLW)  in  the  Mystic  River,  Inner  Confluence,  and  Reserved 
Channel  and  to  a  minimum  depth  of  -40  feet  MLW  in  the  Chelsea 
River  (see  Figure  1) .  The  results  will  supplement  soil  borings 
previously  taken  by  providing  continuous  profile  line  coverage 
of  the  entire  length  of  the  project  area.  The  soil 
classification  analysis  will  enable  accurate  estimates  of 
volumes  of  clay  and  gravel  versus  silt  to  be  made.  In  addition, 
the  extent  of  rock  zones  and  the  volume  or  quantity  of  rock  to 
be  removed  is  desired.  This  will  facilitate^  the  accurate 
positioning  of  any  additional  borings  as  may  be  required, 
particularly  in  areas  of  suspected  rock  outcroppings.  This 
information  will  be  used  to  compute  project  costs  and  in  the 
alternative  disposal  analysis  to  be  included  in  the  project 
Environmental  Impact  Statement. 

Scope  of  Work 

3.  An  array  of  geophysical  instruments  will  be  utilized  to  meet 
the  objectives  of  the  investigation.  The  basic  suite  of  systems 
include: 

a)  a  3.5  to  7.0  kHz  high  resolution  ”pinger”  system, 
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b)  an  integrated,  high  definition  400  Hz  to  5.0  kHz 
"boomer”  system, 

c)  A  new  2.0  to  10.0  kHz  chirper  to  be  used  in  conjunction 
with  the  above  if  it  appears  advantageous  to  the  project,  and 

d)  a  dual  frequency  side  scan  sonar  system. 

This  equipment  will  be  used  to  produce  the  necessary  seismic 
energy  to  obtain  reflection  signatures  from  the  bottom  and 
subbottom  marine  sediments.  Primary  channels  of  acoustic  data 
will  be  acquired  with  a  digital  data  acquisition  system.  Data 
collected  during  the  survey  allows  continual  data  <^ality 
control  and  initial  analysis.  Further  data  analysis, 
processing,  and  visualization  will  be  performed  post-survey  at 
WES.  Acquisition  and  interpretation  of  traditional 
”shades-of-gray”  analog  reflection  records  will  be  performed 
concurrently  with  the  digital  acoustic  impedance  techniques. 

4.  A  dual  frequency  side  scan  sonar  is  typically  towed  during  a 
survey  to  provide  increased  bottom  coverage  and  detect  any 
possible  dredging  or  navigation  hazards.  Although  entirely 
qualitative,  the  side  scan  sonar  can  be  interpreted  to 
differentiate  bottom  sediments  in  terms  of  material  types,  such 
as  sands  versus  clays,  in  a  rapid  continuous  meinner.  In 
addition,  bulkhead  foundation  assessment  and  depths  along  the 
tributaries  can  be  'determined. 

5.  Continuous  subbottom  profiles  of  the  acoustic  reflection 
amplitudes  obtained  from  the  seismic  systems,  an  example  of 
which  is  presented  in  Figure  2,  will  be  provided  within  three 
weeks  after  the  field  survey.  The  analog  reflection  records 
will  also  be  provided  upon  request.  These  records  will  be 
annotated  with  horizontal  control,  existing  core  locations,  and 
preliminary  interpretations  of  the  site  characteristics.  Upon 
receipt  of  ground  truth  information  from  existing  core 
locations,  acoustic  impedance  calculations,  and  density 
estimates  will  be  determined  (see  Figure  3)  to  develop  a  set  of 
seismic  parameters  to  process  the  complete  data  set.  During 
processing,  tabular  listings  of  density  and  material  type  versus 
depth  (see  Figure  4)  as  well  as  plots  of  density  versus  depth 
(see  Figure  5)  will  be  provided  at  sufficient  intervals  along 
the  profile  lines.  The  processed  data  will  be  adjusted  to  MLW 
and  correlated  to  the  positioning  data. 

6.  The  virtually  continuous  data  coverage  of  the  svibbottom 
material  are  incorporated  into  a  terrain  modeling  software 
package  to  provide  a  three-dimensional  (3-D)  view  of  the  project 
area.  An  example  of  3-D  view  of  the  distribution  of  sediment 
densities  in  a  portion  of  the  Gulfport  Ship  Channel,  MS  is  shown 
in  Figure  6,.  Displays  of  this  type  provide  improved  data 
interpretation  and  visualization  assistance  as  compared  to 
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traditional  two-dimensional  illustrations.  Views,  color-coded 
according  to  density  or  material  type,  can  be  created  by 
picturing  the  region  from  any  perspective.  Three-dimensional 
visualization  also  permits  volumetric  calculation  of  the 
(jiffering  sediment  units  and/or  the  amount  of  material  to  be 
removed  through  dredging.  Furthermore,  a  detailed  data  base  is 
established  for  project  monitoring  and  long-term  planning  of 
future  work. 

7.  Proposed  Survey:  Maps  of  the  project  area  with  the  proposed 
survey  lines  illustrated  are  presented  in  Figures  7-10.  The 
survey  is  divided  into  four  parts,  and  is  outlined  as  follows; 

a.  Survey  of  Reserved  Channel;  Improvement  dredging  is 
likely  within  the  existing  Federal  channel  as  well  as  two 
additional  areas  near  the  mouth  of  the  tributary.  Although 
current  proposed  plans  call  for  improvement  to  the  channel  for  a 
distance  of  only  4,000',  survey  lines  will  extend  approximately 
4,500'.  Four  profile  lines  will  be  surveyed,  ^o  lines  will  be 
100'  either  side  of  the  centerline  and  the  remaining  two  will  be 
250'  to  either  side  of  the  centerline.  Four  additional  lines 
with  spacing  of  150*  and  length  of  1,600'  will  be  run  at  the 
mouth  parallel  to  the  the  centerline  starting  400'  to  the 
north.  Finally,  seven  profile  lines  parallel  to  the  centerline 
of  the  35'  main  ship  channel  with  150'  spacing  will  be  taken 
starting  at  200'  downstream  of  buoy  #10  and  proceeding  upstream 
for  3,000'.  Side  scan  sonar  will  be  used  particularly  along 
channel  edges  where  bulkheads  exist.  Specific  effort  will  be 
made  to  identify,  locate,  and  establish  the  approximate  depth  of 
a  power  line  in  the  channel.  Additional  survey  lines  will  be 
performed  depending  on  actual  subbottom  conditions  or  as 
requested  by  New  England  Division  personnel. 

b.  Survey  of  Mvstlc  River;  Improvement  dredging  is  likely 
only  within  portions  of  the  existing  Federal  channel.  However, 
the  entire  channel  between  the  Mystic  River  (Tobin)  Bridge  and 
the  Alford  Street  Bridge  will  be  surveyed.  Seven  profile  lines, 
six  5,500'  in  length,  and  the  northern  most  4,500'  in  length 
will  be  oriented  about  the  centerline  with  spacing  of  200' 
except  for  the  outermost  lines  which  will  have  spacing  of  150 | . 
Additionally,  in  areas  where  rock  is  expected  off  rte  Exxon  Pier 
two' additional  profile  lines  of  1,000'  in  length  will  be  made. 
Side  scan  sonar  will  be  used  particularly  along  channel  edges 
where  bulkheads  exist.  Additional  survey  lines  will  be 
perfomned  depending  on  actual  subbottom  conditions  or  as 
requested  by  New  England  Division  personnel. 

c.  Survey  of  Inner  Confluence;  Improvement  dredging  is 
likely  within  the  existing  Federal  channel  as  well  as  along  the 
northeast  corner  of  the  area.  Because  of  suspected  rock,  eight 
profile  lines,  parallel  with  the  centerline  into  the  Chelsea 
P3_Yer,  will  be  spaced  100'  apart.  Five  lines  will  be  2,200'  and 
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the  other  three  will  be  1,500*  in  length.  Five  lines,  each 
2,000*  in  length,  with  a  spacing  of  150*  will  be  parallel  with 
the  centerline  of  the  Mystic  River  approach.  Finally  two 
additional  profile  lines  will  provide  information  along  the 
northern  and  eastern  edges  of  the  area.  The  northern  edge 
profile  will  be  1,500*  in  length  and  the  eastern  profile  will  be 
3,000*.  Side  scan  sonar  will  be  used  particularly  along  channel 
edges  where  bulkheads  exist.  Specific  effort  will  be  made"  to 
identify,  locate,  and  establish  the  approximate  depth  of  a 
telephone  cable  crossing  the  Inner  Confluence.  Additional 
survey  lines  will  performed  depending  on  actual  subbottom 
conditions  or  as  requested  by  New  England  Division  personnel. 

d.  Survey  of  the  Chelsea  River:  Improvement  dredging  is 
likely  within  the  existing  Federal  channel.  Because  of  the 
narrowness  of  the  channel  below  the  Chelsea  Street  Bridge  only 
two  profiles  will  be  taken  —  one  on  either  side  of  the 
channel.  In  the  upper  reaches  a  total  of  thirteen  profiles  will 
be  obtained.  Directly  upstream  of  the  Chelsea  Street  Bridge 
five  profiles  each  3,000*  in  length  will  be  parallel  to  the 
eastern  bank  and  spaced  150*  apart.  Further  upstreeim  three 
profiles  2,500*  in  length  parallel  to  the  Gulf  Oil  dock  shall  be 
spaced  200*  apart.  Finally,  parallel  to  those  profiles  there 
will  be  four  more  profiles  covering  the  turning  basin  each  being 
1,500*  in  length  and  spaced  200*  apart.  The  last  profile  will 
be  along  the  front  of  the  Global-Gibbs  and  Northeast  Oil  docks. 
Side  scan  sonar  will  be  used  particularly  along  channel  edges 
where  bulkheads  exist.  Specific  effort  will  be  made  to 
identify,  locate,  and  establish  the  approximate  depth  of  utility 
crossings  in  the  channel.  Additional  survey  lines  will  be 
performed  depending  on  actual  subbottom  conditions  or  as 
requested  by  New  England  Division  personnel. 

8.  WES  will  provide  all  necessary  geophysical  equipment  and 
recording  devices  to  complete  the  project.  WES  personnel  will 
set  up  and  operate  all  geophysical  equipment  during  sxibbottom 
surveying  exercises  and  coordinate  survey  operational 
procedures.  WES  also  has  the  authority  to  decide  on  a  daily 
basis,  with  input  from  the  captain  of  the  survey  vessel,  if 
conditions  (weather,  sea  state,  white-capping,  etc.)  are 
acceptable  for  subbottom  profiling  tests.  New  England  Division 
project  personnel  are  encouraged  to  actively  participate  in  the 
data  acquisition  and  processing  exercises  to  both  assist  in  and 
gain  an  understanding  of  the  survey  and  analysis  procedures. 

9.  One  eight-hour  day  for  dock-side  mobilization  to  load,  set 
up,  calibrate,  and,  if  necessary,  troubleshoot  the  geophysical 
equipment  will  be  allowed.  In  addition,  one  day  of  project  site 
reconnaissance  will  be  allowed  to  acquire  subbottom  signatures 
in  all  three  areas  where  physical  cores  have  been  obtained. 

This  will  provide  ground  truth  information  to  calibrate  the 
acoustic  predictions. 
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10.  WES  personnel  will  post-process  all  geophysical  field  data 

and  prepare  a  letter  and  final  report  of  the  test  results.  The 
continuous  subbottom  cross-sections  of  the  acoustic  reflection 
amplitudes  will  be  provided  within  three  ^  ^ 

survey.  These  records  will  be  annotated  with  horizontal 
control,  existing  core  locations,  and  preliminary  interpre  a- 
tions  of  the  site  characteristics.  A  letter  report  of  the 
results  will  be  presented  to  the  New  England  Division  personnel 
irf  one-on-oL  visit  before  mid-January  19 W.  All  aspects  of 
the  survey  and  the  results  detemined  will  be  carefully 
explained  so  that  maximum  benefits  will  be  derived  from  the 
investigation.  The  final  report  will  detail  the  metoods  used  to 
acquire,  calibrate,  and  process  the  data  and  present  the 
findings  and  conclusions  in  a  clear,  concise  manner  that  will 
directly  address  the  purpose  of  the  geophysical  su^ey.  The 
final  report  will  be  completed  by  April  1993  and  will  document, 
in  addition  to  soil  classification  densities  and  volumes, 

such  as  utility  crossings,  bullchead  foundation  information,  and 
suggestions  for  additional  borings. 

finpport  Requirements 

11.  The  New  England  Division  is  expected  to  provide  the 
following  support  for  the  successful  completion  of  ttis  tasK. 

The  cost  of  these  items  is  not  reflected  in  the  project  cost 
estimate. 

a  A  survey  vessel  and  crew  able  to  conduct  the  ^  required 
survey.  The  vessel  must  be  at  least  45*  in  length  with  at  least 
eHSare  feet  of  back  deck  areas  and  an  enclosed  area  to  house 
the  graphic  recorders  and  computers  and  protect  them  from  the 
elements- 

b.  Navigational  and  positioning  support  (horizontal 
control)  to  include  pre-  and  post-plots  of  the  su:^ey  area  with 
reference  marks  every  200'.  If  possible,  the  digital  data 
acquisition  system  will  be  interfaced  directly  with  the  vessel  s 
positioning  system.  An  on-board  x-y  plotter  is  recommended  if 
the  capability  exists. 

c.  Continuous  fathometer  soundings  and  water  depths 
corrected  for  MLW.  This  is  required  to  correlate  depth 
measurements  with  the  geophysical  equipment  and  to  allow 
adjustments  for  tidal  fluctuations. 

d  The  equipment  and  means  necessary  to  take  5  to  10 
gravity  cores  at  locations  within  the  survey  site.  These 
locations  will  be  selected  by  WES  personnel  during  the  d^a 
acquisition  phase  of  the  investigation.  The  insitu  inforaation 
(bulk  density,  mean  grain  size,  water  content,  etc.)  will  be 
used  to  construct  acoustic  impedance  and  density  estimates  for 
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the  local  sediment  regime.  The  cores  will  be  analyzed  by  the 
WES  Soils  Testing  Facility  or  at  the  New  England  Division 
Environmental/Materials  Laboratory. 

e.  If  necessary,  assistance  in  loading  and  unloading  the 
geophysical  equipment  on  and  off  the  survey  vessel. 

12.  It  is  recognized  that  during  data  acquisition,  WES 
personnel  typically  work  12  hour  days  to  minimize  the  number  of 
days  on  the  water  and  to  expedite  data  collection.  The  exact 
days  and  times  will  be  finalized  two  weeks  prior  to  the  start  of 
the  survey.  As  with  all  field  operations,  considerable 
flexibility  should  be  allowed  in  the  schedule. 

Reporting  Results 


13.  The  goal  is  to  have  the  results  available  to  New  England 
Division  personnel  as  soon  as  possible  after  the  field  work  is 
completed  since  the  results  will  likely  be  used  in  preparation 
of  an  Environmental  Impact  Statement  requiring  this  input  by 
15  January  1993.  The  annotated  amplitude  cross-sections  and 
preliminary  data  analysis  will  be  made  availcible  within  two 
weeks  of  completion  of  the  field  survey.  The  results  and 
interpretations  will  be  drafted  into  a  letter  report  and 
presented  to  New  England  Division  personnel  during  January 
1993.  A  complete  and  comprehensive  final  report  of  the  data 
acquisition  and  processing  procedures,  final  results, 
interpretations,  and  conclusions  will  be  completed  by  April 
1993. 
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Figure  2.  Portion  of  a  seismic  reflection  record. 
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During  preparation  of  this  report,  fabrication  of  the  waterborne  integrated 
geophysical  system  (IGS)  was  completed.  That  system,  funded  from  sources 
other  than  the  DRP,  now  includes  a  support  vessel  (shown  in  Figures  76,  77, 
and  78)  housing  an  integrated  GPS  navigation  and  data  management  computer 
with  appurtenances  necessary  for  acquiring  and  processing  a  wide  variety  of 
geophysical  data.  The  vessel,  christened  WATERWAYS  EXPLORER,  was 
designed  to  be  highway  transportable  and  capable  of  stable  operation  in  coastal 
regions,  harbors,  and  inland  waterways.  Needless  to  say,  the  IGS  concept 
greatly  benefited  through  research  performed  under  the  DRP.  (Refer  to 
Appendix  A,  Figure  A1  for  the  block  diagram  depicting  equipment  currently 
on  board  and  including  future  plans  for  increased  capability.)  When  the  vessel 
was  designed,  immediate  emphasis  was  placed  upon  the  incorporation  of 
equipment  (including  drop  tube  sampling  for  immediate  ground  truth  verifica¬ 
tion)  capable  of  acquiring  high-resolution  reflection  seismic  information  com¬ 
patible  with  acoustic  impedance  processing  algorithms.  That  capability  is 
supplemented  with  side-scan  sonar  (with  target  enhancement),  multi-frequency 
Fathometers,  and  magnetometers. 

The  WATERWAYS  EXPLORER  has  surveyed  four  sites  to  date  with  others 
scheduled.  Sites  already  surveyed  include  shorelines,  harbors,  and  inland 
lakes.  Its  tri-pontoon,  shaUow-draft  platform  has  provided  the  necessary  stabil¬ 
ity  for  extremely  high-quality  data  acquisition.  With  continued  use,  however, 
improvements  in  design  will  be  incorporated. 
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